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1.0   INTRODUCTION 

Stainless  steel  type  304,  because  of  its  good  mechanical  properties  and 
corrosion  resistance,  is  often  used  as  structural  and  cladding  material  in 
nuclear  reactors  (1,  2).   In  most  reactors  the  surface  of  stainless  steel, 
either  of  cladding  or  of  structural  material,  is  in  contact  with  water  or  heavy 
water  used  as  coolant  or  moderator  (3).   It  is,  therefore,  imperative  to  study 
the  behavior  of  stainless  steel-304  in  these  two  environments  from  the  view- 
point of  surface  reactions  and  corrosion  resistance. 

Since  corrosion  processes  are  most  often  electrochemical  (4) ,  the  electro- 
kinetic  study  on  a  particular  surface  leads  to  a  clear  conception  of  the  mech- 
anism of  corrosion.   Today  it  can  be  considered  as  an  established  fact  that 
the  primary  cause  of  corrosion  is  the  thermodynamic  instability  of  most  metals 
and  alloys  in  atmospheric  conditions  and  in  water  solutions.   The  ultimate  aim 
of  corrosion  studies  is  to  investigate  methods  of  controlling  corrosion  by 
limiting  the  surface  reactions  either  by  changing  the  environmental  condition 
or  by  changing  the  composition  of  the  material. 

Many  researchers  have  studied  such  problems  by  electrochemical  measure- 
ments on  pure  substances  (5,  6,  7).   Work  on  alloys  by  such  methods,  on  the 
other  hand,  is  lagging.   In  the  case  of  stainless  steel,  e.g.,  most  of  the 
studies  have  been  directed  towards  passivity  (8,  9,  10). 

The  purpose  of  this  thesis  was  to  study  the  electrokinetics  of  a  freshly 
electropolished  surface  of  stainless  steel-304  in  water  and  heavy  water  and  to 
find  the  effect  of  deuterium  substitution  on  the  kinetics  and  mechanism  of 
reactions.   This  work  leads  to  the  determination  of  a  value  for  the  electro- 
lytic separation  factor,  S,  defined  as  the  ratio  of  the  concentration  of 
hydrogen  to  deuterium  in  the  gas  phase  to  the  corresponding  ratio  in  the 
solution  phase  in  an  electrolytic  separation  cell.   Such  a  factor  is  of  extreme 


importance  in  heavy  water  preparation  by  the  process  of  electrolytic  separation, 
which  process  has  received  and  is  still  receiving  utmost  care  in  the  production 
of  heavy  water.   From  knowledge  of  the  separation  factor  and  the  standard  heats 
of  activation  of  the  rate  determining  steps  in  water  and  heavy  water,  further 
insight  into  the  reaction  mechanism  was  obtained. 


2.0  THEORETICAL  DEVELOPMENT 


2.1   Electrode  Potential  and  Overvoltage 

In  a  complete  galvanic  cell  consisting  of  an  electrode  of  metal  M  and  a 
Pt  electrode  in  the  same  solution,  there  are  three  potential  differences: 

1)  V  ,  that  at  the  metal-solution  boundary 

2)  V   ,  at  the  Pt-  solution  boundary 

and   3)   X,  the  contact  potential  difference  between  M  and  Pt. 

If AV  is  the  equilibrium  potential  difference  measured  by  means  of  a 

potentiometer  connected  to  wires  of  composition  M,  then 

AV  =  VM  -  V_   +  X  (1) 

e    Me    Pte 

where  subscript  e  denotes  equilibrium  potentials.   Now  X  =  <t>xi  -   *„  where  4.. 

M    Pt        M 

and  <{>   are  the  work  functions  of  the  metals  M  and  Pt  respectively.   Hence 

AV     =  V        -  V  +  <J>      -   <j>  (2) 

e  Me  Pte        VM       vPt  K   ' 

If,  in  the  cell  considered  above,  the  cathode  carries  a  current  density 

i,  a  different  potential  difference  will  arise  given  by 

AV  -  VM  -  Vpt  +  *M  -  *pt  (3) 

where  AV  ■  irreversible  potential  difference  between  the  metal  cathode  and 

the  Pt  electrode 
and   VM  ■  irreversible  potential  difference  at  the  metal  solution  boundary. 

V-.^.  ■  V__   if  the  platinum  electrode  is  used  as  a  reference  hydrogen 

Pt    Pte       r  J        ° 

electrode  which  is  the  case  in  all  overpotential  measurements. 
AV  -  AV  is  denoted  by  r\,    the  overvoltage.  Or 

AV  =  AV  +  n  .  (4) 

Writing  Eq.  (3)  in  terms  of  the  overvoltage  and  equilibrium  potential,  it 
follows  that 


VM  "  AVe  +  "  +   VPt  -  *M  +  *Pt  "  (6> 

In  the  cell  considered  above,  if  the  Pt  electrode  acts  as  a  hydrogen 
electrode,  then  for  a  given  pH  and  temperature  AV  ,  V   ,  <j>   are  constant. 
Therefore,  Eq.  (6)  can  be  written  as 

VM  =  n  ~  *M  +  K'  (7) 

where  K'  =  AV  +  V   +  <j>   .  (8) 

The  relation  between  the  Cathodic  current  density  i  and  the  potential  V..,  for 

J      c  M 

a  kinetic  electrode  process  in  which  the  reverse  reaction  can  be  neglected,  is 

given  (11)  by  an  equation  of  the  form  (see  below) 

i  =  kn  exp(-aV  F/RT)  (9) 

c     1        M 

where  a  and  k..  are  constants,  F  is  faraday,  R  is  gas  constant,  and  T  is 

absolute  temperature. 

Thus  by  substituting  from  Eq.  (6) 

i  =  k]_  exp(-a(n  "  *M  +  K')F/RT)  (10) 


if  i  =  i  ,  when  n  =  0,  then 
c    o 


i   =  k_  exp(a<}.  F/RT)  (11) 

o    2       M 

and  Eq.  (10)  can  be  written  then  in  terms  of  i  as 

i  =  i  exp(-anF/RT)  (12) 

c    o 

i  ,  known  as  the  exchange  current,  is  characteristic  of  the  particular 

electrode  used  and  is  taken  as  a  basis  for  comparison  of  various  electrodes  in 

overvoltage  studies.   Eq.  (12)  can  also  be  derived  from  the  theory  of  absolute 

reaction  rates  as  follows. 

In  order  for  any  chemical  change  to  take  place,  it  is  necessary  for  the 

atoms  or  molecules  involved  to  come  together  to  form  an  activated  complex. 


This  is  illustrated  by  the  equation 

a  +  bJx*Jc  +  d  (13) 

where  X  is  the  activated  complex.   This  complex  is  regarded  as  being  situated 
at  the  top  of  an  energy  barrier  lying  between  the  initial  and  final  stages, 
and  the  rate  of  the  reaction  is  controlled  by  the  rate  with  which  the  complex 
travels  over  the  top  of  the  barrier. 

According  to  the  theory  of  absolute  reaction  rates  (12)  the  specific  rate 
k1  of  any  reaction  is  given  by  the  expression 

k'  =  T(kT/h)exp(-AF/RT)  (14) 

where  t  is  a  constant  introduced  to  allow  for  the  possibility  that  not  every 
activated  complex  reaching  the  top  of  the  potential  energy  barrier  is  converted 
into  reaction  products,  and  is  called  the  transmission  coefficient.   (Except 
for  limited  cases,  t  is  usually  taken  as  unity.) 

For  a  reversible  electrode  the  specific  rate  of  the  discharge  process, 
i.e.,  of  the  direct  reaction,  may,  therefore,  be  written  as 

k'  =  (kT/h)exp(-AG*/RT)  (15) 

while  for  the  reverse  reaction 

kj;  =  (kT/h)exp(-AG*/RT)  (16) 

the  equilibrium  constant  K  =  k]/k'  =  exp(-(AG  -  AGJ/RT)  (17) 

where  AG  =  standard  fre  energy  of  activation  for  the  direct  reaction 
and   AG_  =  standard  free  energy  of  activation  for  the  reverse  reaction. 
At  the  reversible  potential,  where  Eqs.  (15)  and  (16)  hold  good,  no  current 
passes,  but  if  an  overvoltage,  n>  is  applied,  the  rate  of  the  direct  reaction 
exceeds  that  of  the  reverse  process  and  current  flows  at  a  definite  rate. 
Without  making  any  assumptions  as  to  the  mechanism  of  the  discharge  process, 
it  may  be  supposed  that  the  effect  of  the  excess  potential  is  to  diminish  the 


free  energy  increase  requisite  for  the  formation  of  the  activated  state  from 
reactants,  thus  assisting  the  direct  process,  whereas  the  reverse  process  is 
retarded.   If  the  additional  potential,  n>  acts  across  the  energy  barrier  be- 
tween initial  and  final  states,  and  a  is  the  fraction  of  this  electrical 
potential  operative  between  the  initial  and  activated  states  (a  has  a  most 
probable  value  of  0.5  for  the  discharge  reaction  as  determined  experimentally 
(12)),  then  the  potential,  an,  will  facilitate  the  formation  of  the  latter  by 
decreasing  the  free  energy  by  an  amount  anF.   (See  figure  1) 
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Fig.  1.  Energy  Barrier  Diagram 


At  the  same  time  the  potential  will  oppose  the  reverse  reaction,  increasing 
the  free  energy  required  by  (1  -  a)nF.   The  specific  rates  for  direct  and 
reverse  reactions  at  the  overvoltage,  n,  are  then 


and 


or 


and 


k'  =  (kT/h)exp(-AG  /RT)exp(-anF/RT) 


k£  =  (kT/h)exp(-AG2/RT)exp((l  -  a)nF/RT) 


k'  =  A  exp(-anF/RT) 


k^  =  A2  exp((l  -  a)nF/RT) 


(18) 

(19) 
(20) 
(21) 


where  A  and  A~  are  constants.   Since  no  assumptions  have  been  made  concerning 
the  nature  of  the  reactants  and  resultants,  their  concentrations  may  be  repre- 
sented by  C.  and  C„  units,  i.e.,  atoms,  molecules  or  ions,  per  square  cm  of 
the  electrode  surface  so  that  the  direct  and  reverse  reaction  rates,  in  units 
per  square  cm,  related  to  the  specific  reaction  rate  constants 

vi  -  kici 

v2  =  k;C2 


and 


v  =  AC   exp(-anF/RT) 


V2  =  A2C2  exp((1  "  °0nF/RT)  . 


(22) 
(23) 


The  current  passing  is  determined  by  the  difference  between  these  rates,  and 
if  each  reacting  unit  may  be  regarded  as  carrying  the  equivalent  of  a  single 
charge,  the  current  density  i,  in  amp.  per  sq.  cm,  is  given  by 

i  =  (vx   -   v2)F/N 

=  (C1A1  exp(-anF/RT)  -  Cfa   exp((l  -  a)nF/RT)F/N   (24) 

where  F  is  the  faraday,  i.e.,  96,500  coulombs  and  N  is  Avogadro's  number.  If 
n  is  small,  i.e.,  for  very  low  overvoltages,  the  exponentials  may  be  expanded 
and  all  terms  other  than  the  linear  ones  neglected.   Thus 


ic  =  [CiAi(1  *  anF/RT)  -  C2A2(1  +  (1  -  a)nF/RT)]F/N  .         (25) 

It  is  seen  from  Eq .  (25)  that  at  the  reversible  potential,  where  n  =  0  and 

i   =0,  CnA.,  =  C„A„.   If  this  equality  may  be  assumed  to  hold  also  for  small 
c       1  1    2  2  n     J        J 

values  of  n,  Eq .  (25)  becomes 

i   =  (C1A1nF/RT)F/N  .  (26) 

c     11 

At  low  overvoltages,  therefore,  there  should  be  a  linear  relationship  between 

current  and  overvoltage.   This  has  been  found  to  be  true  for  the  evolution  of 

hydrogen  and  the  deposition  of  metals  (12).   For  higher  values  of  n,  the  rate 

of  the  reverse  reaction  becomes  negligibly  small  in  comparison  with  that  of  the 

discharge  process  so  that  it  is  possible  to  write 

i   =  (C.A.  exp(-anF/RT))F/N  =  i   exp(-anF/RT)  .  (27) 

ell  o 

Current,  i  ,  the  significance  of  which  has  been  indicated  above  is  equal  to 
o 

C  A  F/N.   The  factor,  a,  in  Eq .  (27)  can  be  considered  as  the  product  of  two 

factors,  3  and  A  (11),  where  A  equals  the  number  of  electrons  necessary  so 

that  one  act  of  the  rate  determining  step  can  occur  and  3  is  the  symmetry 

factor.   Eq.  (27)  can,  therefore,  be  written  in  the  form 

i  =  i   exp(-6AnF/RT)  (28) 

c    o 

or  n  =  (RT/3AF)ln  i  -  (RT/3AF)ln  i   .  (29) 

o     -  c 

Eq.  (29)  can  be  written  in  the  form 

n  =  a  -  b  In  i  (30) 

c 

where        a  =  (RT/3AF)ln  i  (31) 

and  b  =  RT/3AF  =  -dn/dln  i   .  (32) 

c 

Eq.  (30)  is  usually  called  Tafel  equation  and  "b"  called  Tafel  slope.   By 

definition,  3  and  A  values  are  given  by 

o<3<1   and    A  =  1  or  2.  (33) 


2.2   Kinetics  of  the  Hydrogen  Evolution  Reaction 

It  has  been  shown  (11)  that  the  possible  reaction  paths  of  the  hydrogen 
evolution  reaction  in  acid  solutions  are  mainly 

H  0+  +  e  +  M   Si°W   MH  ,   +  Hn0  I 

3  ads     2 


MH  ,   +  MH  ,   F^St   2M  +  H0  II 

ads     ads  2 


(A) 


+  Fast 

.  J  +  e  +  M   _>    ....  .  „0, 

3  ads    2 


H„0  +  e  +  M   "    MH  ,   +  H„0        III 


(B) 


(C) 


MH  ,   +  MH  ,   Si°W   2M  +  H0  IV 

ads     ads  I 

Ho0+  +  e  +  M   Si°W  MH  .   +  HLO  V 

3  ads     2 

Ho0+  +  MH    +  e  F_tSt  H0  +  Il.O  VI 

3       ads  2     2 

Ho0+  +  e  +  M   I'^St  MH    +  H„0         VII 
3  ads     2 

and 

Ho0  +  MH    +  e  -°W  H„  +  Ho0        VIII 
3       ads  2    2 

Additional  possible  reaction  courses  are  given  in  reference  (11). 

The  slow  reactions  I,  IV  and  VIII  are  rate  determining  steps  and  are 
commonly  termed  simple  discharge,  atomic  desorption  and  electrochemical  desorp- 
tion  mechanisms  of  the  electrolytic  evolution  of  hydrogen,  respectively. 

Each  of  the  reaction  paths  A,  B,  C,  and  D  is  discussed  below: 
(I)   Reaction  Paths  A  and  B 

The  source  of  the  protons  is  H„0  and  the  atomic  hydrogen  reaction  is 
desorptive. 
General  Kinetic  Equations 

The  reactions  and  velocities  v1  and  v„  to  be  considered  are 

+         Vl 
Ho0  +  e  +  M  +  MH    +  Ho0  (34) 

3  ads    2 
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V2     + 
MH  ,   +  Ho0  +        H„0  +  e  +  M  (35) 

ads    I  3 

V3 
and  MH    +  MH     -*T   H  +  2M  (36) 

ads     ads       2 

The  velocity  of  the  reverse  desorption  reaction  of  Eq .  (36)  is  assumed  to  be 

negligible  at  potentials  greater  than  the  reversible  hydrogen  potential. 

The  potential  difference  between  the  cathode  and  the  solution  represented 
by  AV  consists  of  two  parts.   These  two  parts  are  AVp,  between  the  cathode  and 
the  plane  passing  through  the  center  of  ions  adjacent  to  the  cathode  surface 
and  £,  between  this  plane  and  the  bulk  of  the  solution.   This  is  based  on  the 
concept  of  the  existence  of  an  electrical  double  layer,  at  the  electrode  solu- 
tion boundary  (13),  consisting  of  two  oppositely  charged  layers  of  ions  at  a 
fixed  distance  apart.   This  double  layer  is  regarded  as  equivalent  to  an  elec- 
trical condenser  of  constant  capacity  with  parallel  plates  separated  by  a 
distance  of  the  order  of  a  molecular  diameter. 

From  the  theory  of  absolute  reaction  rates  discussed  above,  the  specific 
rate  constant  of  reaction  (34)  is  given  by 

k'  =  (kT/h)exp[-(AF  +  aAVpF)/RT]  (37) 

where  AF  is  the  free  energy  change  of  reaction  (34)  when  AVp  =  0,  but  since 

the  reaction  rate  v..  ■  a.k'  where  the  activity  of  the  hydrogen  ions  on  the 

electrode  surface  a,  ■  a„+(l  -  x) »  a„+  being  the  activity  of  hydrogen  ions  in 

In  n 

the  electrical  double  layer  and  x  being  the  fraction  of  electrode  surface 
covered  with  adsorbed  hydrogen.  Therefore,  from  Eq.  (37) 

Vl  "  klaH+(1  "  X)e*p(-AVpF/2RT)  .  (38) 

15  2 

If  there  are  10  "  free  spaces  for  adsorption  per  cm  of  the  electrode 

surface  and  x  *-s  the  fraction  of  the  surface  covered  with  adsorbed  hydrogen 

2       15 
atoms,  then  the  number  of  atoms  of  hydrogen  adsorbed,  per  cm  =  x-10   •   Since 

-24 
1  atom  of  hydrogen  is  approximately  10    gm  atom,  then  activity  a„  of  the 
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2      -9. 
adsorbed  hydrogen  in  gm  atoms  per  cm  =  x  10  '    Hence  for  reaction  (35) ,  the 


rate  v~  can  be  expressed  as 


v2  =  k210"9X  exp(AVPF/2RT)  (39) 


and  for  reaction  (36) 


v3  =  k310  18X2  .  (40) 


In  the  steady  state 


Vl  '  V2  "  2v3  =  °  '  (41) 

Solving  Eqs.  (38),  (39),  (40)  and  (41)  for  x  gives  the  real  solution 

-(Ax  +  A2)  +  N/(A1  +  A2)2  +  8AXA3 


4A3 


(42) 


where  A.  =  knau+  exp(-AVpF/2RT)  (43) 

1     1  n 

A2  =  10~9k2  exp(AVpF/2RT)  (44) 

_1  Q 

and  A3  =  10   k3  .  (45) 

For  reactions  (34),  (35)  and  (36),  the  net  current  at  steady  state  is, 

therefore 

i  =  XF(v.  -  v0)  =  2XFv.  =  2FvQ  (46) 

c       1    2        3      3 

and  from  Eqs.  (40),  (42)  and  (46) 

—  1  ft    9         9 
ic  =  2F10  L\3x     =  2FA3xZ  (47) 

=  (F/8A3)[-(A1  +  A2)  +  ((Ax  +  A2)2  +  SA^)1^]2  .      (48) 

(1-1)   Simple  Discharge  is  Rate  Determining  (Path  A) 
(i)   General  condition 

(V1  +  V2)<<V3  (49) 
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since 

VAi 

and 

VA2 

therefore, 

(A1  +  A2)«A3 

or 

10 (A1   +  A2)<A3 

(50) 

where  factor  10  is  selected  as  an  arbitrary  limit  of  significance  (11 )  .   The 
linking  of  v..  and  v„  or  A.  and  A_  by  a  positive  sign  in  the  condition  is  clear, 
since,  by  the  nature  of  the  reactions  concerned,  a  decrease  in  A~  increases  x> 
thereby  decreasing  A.,  and  increasing  A~;  that  is,  a  decrease  in  A„  increases 
the  probability  of  the  slow  discharge  mechanism. 

(ii)  Coverage  of  surface 

The  velocity  v  is  much  greater  than  v„  to  allow  the  slow  discharge 
reaction  to  proceed  in  the  forward  direction,  then 

AX>>A2     or     AX>9A2  .  (51) 

The  factor  9  is  selected  as  an  arbitrary  limit.  Using  conditions  (51)  and  (50) 
in  (42),  it  follows  that 

X  -  (AX/2A3)1/2  .  (52) 

(iii)  Tafel  line 

Substituting  value  of  x  from  Eq.  (52)  and  A1  from  Eq.  (A3)  into  Eq.  (47) 

it  follows  that 

i  -  FA1  (53) 

c     1 

or  i  -  FK^^  exp(-AVpF/2RT)  .  (54) 

c      ill 

Differentiating  AVp  of  Eq.  (54)  with  respect  to  In  ic>  it  follows  that 

-dAVp/d  In  i  -  2RT/F  .  (55) 

Also  it  can  be  shown  that 

b  -  -dVp/d  In  i  -  -dn/d  In  i  -  RT/3XF  .  (56) 

c  c 
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Therefore,  when  simple  discharge  is  a  rate  determining  step 

b  =  2RT/F 

a  -  X0  -  1/2  . 
Since  A  =  1,  0  =  1/2  .  (57) 

(1-2)   Atomic  Desorption  is  Rate  Determining  (Path  B) 
(i)   General  condition 
Comparing  with  (50),  the  condition  is 

10A;}<(A1  +  A2)  .  (58) 

(ii)   Coverage  of  surface 

Following  the  same  idea  used  in  obtaining  Eq.  (51) 

9A2<AX  .  (59) 

From  Eq.  (42)  the  value  of  x  Is 


-(A1  +  A  )  +  N/(A1  +  A2)2  +  8AXA3 

-4A^ (42) 


[-1  +  (1  +  (8A1A3)/(A1  +  A2)2)1/2] 
(AA3)/(A1  +  A2) 


From  condition   (51),   (8A  AJ/(A     +  A2)<1,   hence 


X  -  -(A1  +  A2)/(4A3)   +   (A1  •+  A2)(l  +   WA^)/ 

(Ax  +  A2)2)/4A3 


or 


(60) 


(61) 


X  =  A1/(A1  +  A2)         ,  (62) 

(iii)      Tafel   line 

From  Eqs.    (43),    (44),    (47)    and    (62) 

ic  -  2F10"18k3(A1/(A1  +  A2))2  (63) 


u 


=   2Fk310        [1/(1  +    ((k210   y)/k1aH+)exp(AVpF/RT)r  (64) 

—  1  ft  — Q  — ? 

=  2Fk310  x°[l  +  ((k210   )/k1aH+)exp(AVPF/RT)]  L    .      (65) 

-18  -9 

Therefore,     In  i   =  ln(2Fk„10   )-2  ln[l  +  ((ko10   )/k.au+)exp(AVpF/RT) . (66) 

C  J  /         1  H 

Evaluation  of  d(AVp)/d  In  i   yields 

1  +  [(k210~9)/(k1aH+)]exp(AVpF/RT) 
-d(AVp)/d  In  ic   =  [ (2Fk010-y ) / (kn  RTa+) ] exp (AVpF/RT)  '  (67) 

From  Eq.  (58),  it  is  seen  that  the  values  of  a  and  3  are  complex  functions  of 
potential  and  can  be  calculated  from  (58)  by  the  use  of  (23)  and  (24),  taking 
X  as  2  (since  two  electrons  are  involved  for  Path  B) . 

Two  important  limiting  conditions  arise  from  Eq.  (67). 
Limiting  condition  1 

If  [(k210~9)/(k1aH+)]exp(AVPF/RT)>10  .  (68) 

Then  from  Eq.  (6  7) 

-d(AVp)/d  In  i  =  RT/2F  =  b  (69) 

where  b  is  the  Tafel  slope.   Therefore,  comparing  with  the  general  form  of 
Tafel  slope,  b  =  RT/A3F,  a  is  2  and  3  is  1  because  A  =  2.  (70) 

Limiting  condition  2 

If  AVp  tends  to  -°°,  b->-°°  (71) 

ct->o  and  3+0  .  (72) 

Inclusion  of  the  limiting  condition  (68)  in  (65)  gives 

i   =  2Fk„[(k.au+)/kJ2exp(-2AVpF/RT)  (73) 

C        J    1  n      Z 

and  of  (71)  in  the  same  equation,  gives 

—  1  ft 

i  =  2Fk.,10    .  (74) 

c      J 
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II   Reaction  Paths  C  and  D 

The  source  of  protons  is  HO  and  the  electrochemical  reaction  is  desorp- 
tive. 
General  kinetic  equations 

The  reactions  and  velocities  to  be  considered  are 

+         Vl 
H_0  +  e  +  M  -►   MH    +  Ho0  (75) 

3  ads    2 

V2     + 
MH    +  Ho0   -v   H„0  +  e  +  M  (76) 

ads    2         3 

+  V4 

and  Ho0  +  MH  ,   +e  -*   H„  +  M  +  Ho0  .  (77) 

3       ads  2.2 

v..  and  v„  are  given  by  (38)  and  (39).   v,  is  given  by 

v4  =  k410~9aH+x  exp(-AVpF/2RT).  (78) 


At  the  steady  state 


Vl  "  V2  "  V4  =  °  (79) 

From  Eqs.  (38),  (39)  and  (78)  value  of  x  is 

X  -  A1/(A1  +  A2  +  A4)  (80) 

where  A  =  k,10"9a  +  exp(-AVpF/2RT)  .  (81) 

The  current  is    i    AFv, 

c      4 

=  AFk,10~9au+x  exp(-AVPF/2RT)  (82) 

H  n 

=  2Fk.l0~9au+x  exp(-AVpF/2RT)  .  (83) 

From  Eqs.  (43),  (44),  (80)  and  (83) 

-9 
2Fk,10   au+  exp(-AVpF/2RT)k.au+  exp(-AVpF/2RT) 

h H 1  H 

(k.au+  +  k.l0~  au+)exp(-AVPF/2RT)  +  10~9k„  exp(AVpF/2RT) 
1  n      h  rl  i. 
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2Fk.kA(a  +)210   '    exp(-AVPF/2RT) 


(klaH+  +  k41U   9aH+)    +   i0  9|c2   exP(AVPF/kT> 


and  from  Eq .  (84) 


(84) 


d(AVp)     1  +  [(k210~9)/(k1  +  kA10  9)aH+]exp(AVpF/RT) 

d  ln  Xc    -F/2RT[1  +  [(k.,10_9)/(k1  +  k.lO~9)au+]exp(AVpF/RT)] 

2  1     h  n 

-L(Fk  10"9/(k  +  k  l()"9)a  +RT)exp(AVpF/RT)J  .  (85) 


(II-l)   Simple  Discharge  is  Rare  Determining  (Path  C) 
(i)   General  condition 

:()(A,  f  A,  )<A.  (86) 

(ii)   Coverage  of  surface 
From  Eqs.  (80)  and  (86) 

X  -  Aj/A^  . 

From  Eqs.  (43)  and  (81) 


(87) 


X  -  (k^/UO-9^)  .  (88) 

(iii)  Tafel  line 

Condition  (86)  when  put  in  Eq .  (84)  gives 

i  »  2FkiaH+  exp(-AVpF/2RT)  .  (89) 

This  relation  is  very  similar  to  tiiat  found  foi  a  rate  determining 
discharge  reaction  followed  by  an  atomic  hydrogen  deaotption  step,   in  Eq,  (89) 

0.S  (*U) 

and  .-.lin A:  ■' 

;■    1/4  .  v:M) 

(Ii-2)   Kinetics  When  Electrochemical  Step  is  Rate  Determining  (Path  D) 
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(i)   General  condition 

10A4<A1  +  A2  (92) 

(ii)   Coverage  of  surface 

Special  case  (a)  10A  <A  (93) 


From  Eqs.  (80),  (92)  and  (93) 


X  -  A1/A2  (94) 


Special  case  (b)  10A  <A  (95) 

From  Eqs.  (80),  (92)  and  (95) 

X  -  1  (96) 

(iii)   Tafel  line 
Special  case  (a).   Conditions  (92)  and  (93)  and  Eq .  (84)  give 

i  =  2F(k1/k2)k4(aH+)2exp(-3AVpF/2RT)  (97) 

therefore,  a  =  3/2 

and        A  =  2;  6  =  3/4   .  (98) 

Special  case  (b) .   From  conditions  (92)  and  (95)  and  Eq.  (84),  it  follows 

i  =  2Fk,a  +10   exp(-AVpF/2RT) .  (99) 

Therefore        a  =  0.5;  0  =  1/4  as  X  =  2  .  (100) 

Results  of  the  derivations  presented  above  for  various  reaction  mechanisms 
are  given  in  Table  1. 


Table  I.   Characteristics  of  various  mechanisms  of 
cathodic  hydrogen  evolution  in  acid  solutions 
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With 

RT 
b  "  AgF 

Mechanism 

Conditions  for 
Application 

Symmetrical 
Energy  Barrier 

A 

a 

B 

Slow  Discharge 

10 (Ax  +  A2)<A3 

Fast  Atomic  Desorption  (A) 

9A2<A1 

1 

1/2 

1/2 

2RT/F 

Slow  Atomic  Desorption 

10A.<A  +  A 

Fast  Discharge  (B) 

9A2<A1 

2 

2 

1 

RT/2F 

Slow  Discharge 

Fast  Electrochemical  (C) 

10 (A1  +  A2)<AA 

2 

1/2 

1/4 

2RT/F 

Slow  Electrochemical 

A,  +  A„>10A. 
12     4 

Fast  Discharge  (D) 

A2>10A1 

2 

3/2 

3/4 

2RT/3F 

A1>10A2 

2 

1/2 

1/4 

2RT/F 
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2.3   Stoichiometric  Number,  p 

The  stoichiometric  number  is  defined  as  the  number  of  acts  of  the  rate 
determining  step  accompanying  one  act  of  the  overall  hydrogen  evolution 
reaction.  . 

The  hydrogen  evolution  reaction  may  be  represented  in  acid  solutions  by 
the  overall  reaction 

2H  0+  +  2e  t   H   . 

If  there  is  one  step  in  the  overall  reaction  which  is  rate  determining,  then 
the  rate  v  of  the  forward  direction  of  this  step,  derived  from  Eq.  (37),  is 

v±   =  a1(kT/h)exp[-(AG*  +  2B(AVc  -  0(F/u))/RT]  (101) 

where  a  =  product  of  the  activities  of  reactants  of  the  rate  determining  step 

at  the  site  of  reaction. 
The  analogous  equation  for  the  reverse  reaction  is 

v2  =  a2(kT/h)exp[-(AG2  -  2(1  -  B)  (AV"c  -  C)(F/u))/RT]  .      (102) 

The  currents  i..  and  i„  corresponding  to  v..  and  v„  are 

i1  =  2ev1/u 

±2  =   2ev2/u  (103) 

where  e  is  the  electronic  charge.   When  AV  =  AV„,  i,  =  i„  =  i  .   In  the 

c     R   1    2    o 

general  case,  since  AV  +  n  ■  AV  and  i  =  i   -  i  ,  it  follows  from  Eqs.  (101), 
(102),  and  (103)  that 

i  =  i  [exp(-2BnF/uRT)  -  exp(2(l  -  B)nF/uRT)]  .  (104) 

At  values  of  n  greater  than  zero,  Eq.  (104)  reduces  to 

i  =  i   exp(-2BnF/yRT)  (105) 

c    o 

and  when  n  -*■  o,  Eq .  (104)  approximates  to 
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i  -  -2i  dF/mRT  .  (106) 

c      o 


Differentiating  and  rearranging  Eq.  (106), 


u  -  (-21  F/RT) (dn/di  )   .  rt  .  (107) 

o         c  n  ■*■  o 
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2.4   Deuterium  Substitution 

Several  theories  on  the  electrolytic  hydrogen  deuterium  separation  from 
aqueous  acid  or  alkaline  solutions  have  attributed  the  difference  in  rates  of 
the  electrolytic  production  of  H  ,  HD  or  D  primarily  to  the  different  zero- 
point  energies  of  the  OH  or  OD  bonds  undergoing  dissociation  in  the  proton  or 
deutron  transfer  in  the  first  simple  discharge  step,  i.e., 

h\0+  +   e  +  M     -*■     MHJ      +  Ho0 
3  ads  2 

in  the  overall  hydrogen  evolution  reaction  at  the  electrode  surface.   This 

procedure  is  an  over  simplification  and  leads  to  values  of  the  separation 

factor,  S,  which  are  not  in  good  agreement  with  experiment  for  various  metals 

and  fails  to  explain  the  clear  division  of  values  of  S  determined  for  various 

metals  in  acid  solutions  into  two  groups. 

In  order  to  explain  some  of  the  differences  between  the  theory  and 

experiment,  one  or  more  of  the  following  factors  are  to  be  considered. 

1.  Differences  of  hydration  energy  of  the  isotopic  ions. 

2.  Relative  concentrations  of  the  discharged  ions. 

3.  Difference  of  ionization  potential  of  D  and  H. 

4.  Differences  of  entropy  of  activation  for  deuterium  and  hydrogen 
evolution. 

5.  The  rate  determining  mechanism  involved  in  the  overall  reaction. 

6.  Whether  H„  and  HD  are  produced  in  the  rate  determining  step  or  not. 
By  considering  all  the  quantities  which  determine  the  heat  contents  of 

the  initial,  activated,  and  final  states  of  various  possible  rate  determining 
steps  in  the  overall  electrolytic  H~  ,  HD  and  D?  evolution  reactions,  it  is 
shown  that  some  important  experimental  facts  can  be  explained. 
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2.4.1  Difference  in  Zero  Point  Energy 

The  major  factor  which  contributes  to  the  free  energy  difference  is  the 
difference  in  zero  point  energy  between  a  bond  to  deuterium  and  the  corres- 
ponding bond  to  hydrogen.   The  potential  energy  curves  (Morse  curves)  for  a 
bond  to  deuterium  and  a  corresponding  bond  to  hydrogen  are  essentially 
identical. 


E£  hydrogen  zero  point  energy 
E0  deuterium  zero- point  energy 


INTERATOMIC         DISTANCE 


Fig.  2.  Morse  Curve  Relating  to  Potential  Energy 
to  Interatomic  Distance 
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The  lowest  energy  level  for  any  bond  corresponds  to  (l/2)hv,  where  h  is 
Planck's  constant.   This  lowest  energy  level  is  known  as  the  zero  point  energy. 
This  corresponds  to  the  vibrational  energy  of  the  bonds  at  absolute  zero.   At 
room  temperature,  approximately  99%  of  the  bonds  are  in  this  vibrational  energy 
level.   There  is  a  difference  in  zero  point  energy  for  a  bond  to  hydrogen  and 
the  corresponding  bond  to  deuterium.   This  arises  from  the  effect  of  the 
difference  in  mass  on  the  stretching  frequencies.   The  shape  of  the  bottom  of 
the  energy  curve  governs  the  force  constant  for  the  stretching  vibration,  and 
this  is  related  to  the  frequency  by  Hooke's  law 

v  =  (1/2tt)(K  /p  )1/2  (108) 

F   m 

where  K„  is  the  force  constant  and  u   is  the  reduced  mass,  which  is  approxi- 

F  m 

mately  equal  to  1  and  2  for  most  hydrogen-  and  deuterium-containing  bonds 
respectively. 

The  difference  in  zero  point  energy  has  two  consequences.   The  first  of 
these  is  a  difference  in  dissociation  energy  for  the  hydrogen-  and  deuterium- 
containing  bonds.   The  dissociation  energy  is  the  difference  in  energy  between 
E   and  the  zero  point  energy  in  Fig.  2.   Since  the  deuterium  compound  has  lower 
zero  point  energy  than  the  corresponding  hydrogen  compound,  the  deuterium 
compound  will  have  a  larger  dissociation  energy,  and  thus  be  more  stable  than 
its  hydrogen  analog.   The  remaining  consequence  can  be  explained  by  noting  Fig. 
3  in  which  the  potential  energy  of  the  system  is  plotted  against  the  distance 
along  the  reaction  coordinate.   Following  the  idea  of  activated  complexes 
discussed  in  Section  2.1,  if  one  assumes  the  bond  undergoing  reaction  to  be 
relatively  weak  in  the  activated  complex  in  comparison  to  the  bond  in  the 
reactants,  the  effect  of  zero  point  energy  on  the  rate  of  reaction  becomes 
apparent.   The  weak  bond  in  the  activated  complex  reflects  a  low  force  constant, 
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activated 
complex 


products 


REACTION       COORDINATE 


Fig.  3.  Potential  Energy  Profile 


and  since  the  difference  in  zero  point  energy  decreases  with  decreasing  force 
constant,  the  difference  in  zero  point  energy  for  the  bond  in  the  activated 
complex  will  be  small.   Thus,  the  difference  in  zero  point  energy  in  the 
reactants  will  result  in  a  larger  energy  of  activation  for  the  deuterium  con- 
taining compound.   in  this  reasoning  it  is  assumed  that  the  bonds  in  the  mole- 
cule which  do  not  participate  in  the  reaction  are  not  affected  during  the 
process.   This  appears  to  be  a  good  approximation  in  most  cases. 

From  the  consideration  of  free  energy,  it  would  then  appear  that  hydrogen- 
containing  bonds  would  undergo  reaction  with  more  ease  and  at  more  rapid  rates 
than  corresponding  deuterium  bonds. 
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2.4.2   Separation  Factor 
The  separation  factor,  S,  is  defined  by 

s  =  (c  V)  (109) 

where  C  =  concentration  of  hydrogen 

Cn  =  concentration  of  deuterium 

II  =  after  the  separation 
and    I  =  before  the  separation. 

In  electrolytic  separation,  II  refers  to  the  gas  phase  and  I  refers  to 
the  solution  phase. 

The  ratio  of  H  to  D  in  the  gas  is  related  to  the  ratio  of  the  currents 
for  hydrogen  and  deuterium  ion  neutralization  at  a  given  metal-solution 
potential  difference. 

During  the  electrode  process  of  separation  of  H  and  D,  there  is  a  definite 
change  in  free  energy  which  can  be  written  (14)  as 

AF  =  AF°  +  RT  In  K  (110) 

where  AF  =  change  in  free  energy  of  the  reaction 

AF°  =  change  in  free  energy  of  the  standard  state  of  unit  activity 
R  =  gas  constant 
T  =  absolute  temperature 
and    K  =  reaction  constant. 
For  the  case  of  equilibrium  AF  =  0;  therefore,  -AF°  =  RT  In  K  or 

K  =  exp(-AF°/RT)  .  (Ill) 


tt  (Products) 

-~~ — — ,=  exp(-AF°/RT)  . 

Tf  (Reactants)     r 


Therefore,   K  =  -~- — —  =  exp(-AF°/RT)  .  (112) 


a 


For  the  hydrogen  reaction,  since  tt  (Products)  =  (C„)TT 

a  nil 
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then,  (CJTT    =   7i    (H   reactants)exp  (-AF°  /RT)     .  (113) 

11113  n 

For  the  deuterium  reaction 

(C  )    =  tt  (D  reactants)exp(-AF°/RT)  .  (114) 

U      L  L  EL  U 

In  the  case  of  electrolytic  separation,  the  change  in  the  free  energy  of 
standard  state,  i.e.,  AF° ,  is  due  to  the  change  in  the  free  energy  of  activa- 
tion for  the  rate  determining  H  or  D  producing  reactions  and  is  denoted  by 
(AG°)  .   Dividing  (113)  by  (114)  and  writing 

AF°  =  (AG°)*  (115) 

it  follows  that 

tt  (H  reactants)  ^        ^ 

<WlI  ■  i"(D  reactants)  °»PK<M£>  "  <AGH>  >/RT'  '  (116> 

a 

Derivation  of  the  separation  factor  using  partition  functions. 

The  separation  factor,  S,  which  is  due  to  the  different  energies  of 
molecules  can  also  be  evolved  from  the  point  of  view  of  partition  functions. 
The  partition  function,  Q,  is  defined  by  the  relation 

Q  =  g  e   o    +  g  e   1    +  g  e   2    +  ...  (11/) 

where  e  ,  e  ,  and  e   are  the  energies  of  a  single  molecule  in  the  various 
possible  levels  (15,  16)  and  g  ,  g  ,  are  constants,  or 

-e./KT 

Q  =  E  g.e   X     .  (118) 

i  &i 

The  entropy,  s,  is  expressed  in  terms  of  the  partition  function  (16)  as 

s  =  E/T  +  R  ln(Q/N)  +  R  (119) 

where  E  =  total  energy 

and    N  =  Avogadro's  number. 

Using  the  thermodynamic  definition 
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F'  =  E  +  PV  -  Ts  (120) 

=  E  +  RT  -  Ts  (121) 
and  combining  Eqs.  (119)  and  (121),  it  follows  that 

F'  =  -RT  ln(Q/N)  .  (122) 

The  partition  function  Q  is  related  to  the  partition  functions  based  on  the 

energy  zero  by  the  relation  (16) 

-E  /RT 
QE  =  Qe  °  (123) 

o 

Hence  Eq .  (122)  can  be  written  in  terms  of  Q„  where  E   is  the  total  energy 

b  O 

o 
based  on  energy  zero 

F  =  -RT  ln(Q/N)  +  E   .  (124) 


Using  the  superscript  "o"  when  the  substances  are  in  their  standard 
states 

F°  =  -RT  ln(Q°/N)  +  E°  .  (125) 

In  a  reaction  aA  +  bB  ~Z   1L  +  mM 

AF°  =  F°        -  F° 

Products    Reactants 

=  -RT[ln(Q°/N)1  +  ln(Q°/N)m  -  ln(Q°/N)a  -  ln(Q0/N)b]  (126) 

+   E°  -  E° 

o  Products    o  Reactants 

or  AF°    -   -RT        In    K'/NhV/N)}}        +   A£o  (12?) 

(Q°/N)a(Q°/N)£ 

where  AE°  =  E°  -  E°  .   Also 

o    o  products    o  reactants 

-AF°  =  RT  In  K  .  (128) 

Equating  (127)  and  (128) 

RT  In  K  =  RT  In    W^L  <QB/N)g   _  AE°  (129 

(Q°/N)^(Q°/N)g       ° 
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K  «    K'/NhV/N)}?   e-AEVRT    _  (13Q) 

(Q°/N)^(Q°/N)g 

* 
In  the  case  of  formation  of  activated  complex,  X  ,  by  the  reactants 

A  +  B  t   X 
(Q°/N)  * 

K  =  WmTJovm:  -p(-AEo/RT)  (131) 

A       B 

where  AE°  =  heat  of  activation  at  absolute  zero  =  AH  because  APV  =  0  and 
o  o 

because  K  =   tt    (Products) /ti    (Reactants).  (132) 

a  a 

The  above  equations  are  now  used  for  H  and  D  reactions.   From  Eq.  (130) 

tt  (Products)  =  (C  ) 
a  H  11 

(Q°/N)x*  ^ 

=  tt  (H  reactants)   r  ,no  /xn 7   exp(-AH  U/RT)  .       (133) 

a  tt L  (Q  /N)u        .  J        o,H 

H  reactants 

Similarly, 

(Q°/N)  *  t 

(VlI  =  VD  '««««.)  ,[(Q./M)  f  exp(-4H    /RT)  .        (134) 

D  reactants 

Therefore, 

tt  (H  reactants)   [F°/tt-(H  reactants)]        ^       ^ 

'Wll  =  /(D  reactants)   [F°/tt  (D  reactants)]  exp[(AHo,D  "  AHo,H)/RT] 
3.  jj   r 

tt  (H  reactants)  [F°/tt  (H  reactants)] 
or      S  =  (C/C„) 


'D  IT  I   tt  (D  reactants)  [F°/tt..(D  reactants)] 
a  D   r 


exp[(AH*jD-  AH*)H)/RT] 


(135) 


where  tt  =  the  product  of  partition  functions  of  H  and  D  containing  reactants 

F°  =  total  partition  function  for  the  activated  state  involving  H 
H 

F°  =  total  partition  function  for  the  activated  state  involving  D 
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2.4.3  Proton  and  Deuteron  Affinities  of  HO  and  DO 

The  proton  affinity,  i.e.,  the  enthalpy  change  in  the  reaction 

H+   +  H„0    ■>  Ho0+ 
gas     2  gas     3  gas 

is  -182  Kcal/mole  (17). 

The  deuteron  affinity  of  H„0  and  D„0  differs  from  the  proton  affinity  of 
water  on  account  of  the  different  zero  point  vibrational  energies  of  the  bonds 
in  DO  and  HO  and  D  0+,  H  D0+,  and  H  0+. 

The  difference  of  proton  and  deuteron  affinities  of  H„0  is  represented  by 
AM/r+   u+>.  (18)  and  is  given  by 

AH(D+   -   H+)    "    [(!(1/2)hVo}DH20+  "    4(1/2)hV°V    ' 

(136) 
6  3 

[(£(l/2)hv   )        +  -    (E(l/2)hv   )         ]    . 
]_  o    HO  2  °    H„U 

The  terms  are  evaluated  from  the  spectroscopic  ground  state  frequencies 
for  the  6  vibrational  modes  of  the  oxonium  ion  and  three  of  H„0  molecule.   The 
number  of  modes  of  vibration  are  given  by  (3n  -  6) ,  where  n  is  the  number  of 
atoms  in  the  ion  or  molecule. 

After  substituting  the  values  for  the  zero  point  energies  reported  by 

Dennison  (19)  Eq.  (136)  gives 

AH/^+   u+n  ■  (19.60  -  12.70)  -  (21.40  -  12.70)         (137) 
v.D  -  H  ; 

=  -1.8  Kcal/mole. 
In  the  presence  of  excess  of  H„0  at  the  start  of  a  normal  separation 
experiment,  both  the  ions  HO  and  DH„0  will  be  almost  exclusively  solvated 
by  H_0  molecules,  and  the  difference  of  overall  solvation  energy  of  D  and  H 
in  excess  H„0  will  be  equal  to  the  difference  of  deuteron  and  proton  affinity 
of  H?0.   Since  HO  and  DH  0  have  the  same  diameter  (the  mean  intermolecular 
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distances  between  0  and  H,  and  0  and  D  atoms  in  the  OH  and  OD  bonds  differ  by 
less  than  0.001  A°  in  the  ground  state)  and  charge;  hence  the  solvation 
energies  of  H^O   and  DH20  will  be  almost  identical. 

2.4.4   Solvation  Energy  of  Ions 

Solvation  energy  of  ions  in  a  medium  of  dielectric  constant  d  may  be 
regarded  as  equivalent  to  the  difference  in  the  electrostatic  energy  of  a 
gaseous  ion  and  that  of  an  ion  in  the  medium.   In  order  to  evaluate  this 
quantity  theoretically,  use  is  made  of  the  method  by  Born  (20);  the  free 

energy  increase  accompanying  the  charging  of  a  single  gaseous  ion  in  medium 

2  2 
of  dielectric  constant  unity,  is  z  e  /2r,  (13),  where  ze  is  the  charge  carried 

by  the  ion  and  r  is  its  effective  radius,  the  ion  being  treated  as  a  conducting 

sphere.   If  the  same  ion  is  charged  in  a  medium  of  dielectric  constant  d,  the 

2  2 
free  energy  change  is  z  e  /2dr,  and  so  the  increase  of  free  energy,  AG, 

accompanying  the  transfer  of  1  gm  mole  of  the  gaseous  ion  to  that  medium, 

which  is  the  free  energy  of  solvation,  is  given  by  the  Born  equation 

AG  =  [(-Nz2e2)/(2r)](l  -  1/d)  .  (138) 

+  + 

The  difference  of  overall  solvation  energy  of  H   in  H90  and  D   in  D^O  can 

also  be  obtained  (21)  from  the  electromotive  force  of  the  cell, 


Pt  (H2) 


HC1  (IN)  in  HO 


DC1  (IN)  in  D20 


D2  (Pt) 


which  is  -0.002  V  (21) . 

To  estimate  the  liquid  junction  potential  in  the  above  cell,  use  is  made 

of  the  Henderson  equation  (13) , 

(U   -  V,)  -  (U   -  V  ) 

ET  =  (RT/F)— ~ —  ln[(U'  +  V')/(u'  +  v'  )]  (139) 

L         (U'  +  V  )  -  (U'  +  V')       I    1     2     Z 
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U1  =  ^(C+U+)1   ;    V±   =  Z(C_U_)1 

Uj  =  2(C+Z+U+)1     Vj  =  E(C_Z_U_)1  (140) 

where  C  and  C_  refer  to  the  concentrations  of  the  cations  and  anions  respec- 
tively, in  gram  ions  per  litre,  U  and  U_  are  the  corresponding  mobilities 

and  Z,  and  Z_  their  valencies  and  where  the  suffixes  1  and  2  refer  to  the  ions 
+ 

in  solutions  1  and  2  respectively. 

In  the  present  case,  solution  1  is  a  HC1-H  0  solution  and  solution  2  is  a 

DC1-D  0  solution,  Z  =  Z_  =  1,  hence 

U  =  u1   :   V  =  V' 
1    1  *        1    1 

U2  =  U'   ;   V2  =  V'   .  (141) 

Eq.  (139)  reduces  to 

(U   -  V  )  -  (U   -  V  ) 
EL  =  (RT/F)  ,(u1  +  /)  _   '  +     '      InKl^  +  Vl)/(U2  +  V,)]     (142) 


Ul  =  W   ;   U2  =  W 


Vl  =  CcrUCl-   ;   V2  =  CC1-UC1-   '  (143) 

Since  the  HC1-H-0  and  DC1-D-0  solutions  are  of  the  same  concentrations, 
it  is  assumed  that 

CH+  =  CC1" 

V  =  CC1- 

CH+  =  CD+   .  (144) 

Eq.  (142)  then  reduces  to 


EL  =  (RT/F)ln[(U1  +  V1)/(U2  +  V^] 


=  (RT/F)ln[(CH+UH+  +  CC1_UC1_)/(CD+UD+  +  CC1_UC1_)] 

=  (RT/F)ln[UR+  +  UC1_)/(UD+  +  U^-)]  .  (145) 
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Defining  transference  numbers  for  hydrogen  and  deuterium  ions  as 

(t+)H  -  V./(uH+  +  uci-5 

(t+)D  =  UD+/(UD+  +  UC1_)  (146)' 

and  substituting  Eq.  (146)  into  Eq.  (145)  it  follows  that 

EL  =  (RT/F)ln[(UH+/(t+)H)/(UD+/(t+)D)]  .  (147) 

Using  the  mobility  data  given  in  references  (22)  and  (23)  and  the  trans- 
ference number  data  given  in  reference  (24) 

E  =  .007  V.  (148) 

The  corrected  E.M.F.,  therefore,  is  -.009  V  and  corresponds  to  a  standard  free 

energy  change,  AG  ,  given  by 

^r  rv        -0.009  x  96500    -  -.  _  .,      .  ,,/cn 

-AG  =  nFE  =  nnnn ; — -,0-,  =   -0.21  Kcal/mole  .         (149) 

o  1000  x  4.183 

This  can  be  equated  approximately  to  the  heat  content  change  in  the  cell 

reaction,  AG°  cell. 

The  difference  of  overall  heat  of  solvation  A[AH     +      +]  of  D   in 

D„0  and  H   in  HO  is  then  given  by 

-i(AHs,D30+  -  H30+>  "  M(1/2)DJ2  .  ^    +  MID  .  „]  +  AG'  cell(150) 

where  A[(1/2)D°      ]  is  the  difference  between  the  bond  dissociation  energies 

D2  '  H2 

of  D»  and  H  ,  A[I     ]  is  the  difference  between  the  ionization  potentials  of 
2.  L  D  -  H 

D  and  H  atoms. 

Using  the  values  shown  in  Table  II  and  the  value  of  AG°cell  calculated 

above  gives 

A[AH  _  n+       „   _+]  -  -1.2  Kcal/mole  .  (151) 

s,D„0   -  H.,0 
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2.4.5   Rate  Determining  Steps 

As  indicated  in  Section  2.2,  the  following  reactions  have  been  suggested 
as  rate  determining  steps  in  hydrogen  evolution  at  various  metals  on  the  basis 
of  modern  experimental  data  obtained  in  ultra-purified  solutions. 

I   H...0  +  M  +  e~  ->  H„0  +  MH  ,    :  Simple  discharge  or  slow  discharge 
J  I  ads 

II   H_0  +  MH  ,,  +  e~  -*  M  +  Hn  +  H„0   :  Electrochemical  desorption  or 
3       ads  2     2 

atom  ion  desorption 

III   MH  .   +  MH  ,   ->-  H„   :  Atom  desorption  or  atom  recombination 
ads     ads    2 

Similarly  analogous  reactions  are  suggested  for  partially  and  completely 

deuterated  species. 

Case  I:   Simple  Discharge 

(a)   Simple  discharge  in  water  containing  the  normal  isotopic  mixture  of 
H  and  D. 

The  isotopic  molecules  are  H„0,  HDO  and  D„0  and  oxonium  ions  are  H„0  , 
D  0+,  H  D0+  and  HD  0+. 

The  D  content  of  normal  water  is  0.015  mole  %  and  the  equilibrium  constant 
for  the  reaction 

HO  +  DO  -*  2H0D 

is  3.8  at  25°C  (25)  which  gives  the  equilibrium  activities  of  the  isotopic 

water  species  as  follows 

a..  ..  =  55  mole/liter 
H2° 

aH0D  =  l   X  10"2 

aD  Q  =  5  x  10~7    .  (152) 
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For  calculation  of  S,  the  separation  factor,  from  Eq .  (135)  the  ratios  of 
pairs  of  equilibrium  activities  of  isotopic  oxonium  ions  are  required.   These 
ratios  can  be  calculated  from  the  following  equations  if  the  ratios  of  the 
respective  equilibrium  constants  can  be  obtained. 

aH30+  =  KiaH+aH20 

aDH20+  =  K2aD+3H0D 

aD2HO+  =  K3aD+3H0D 

aD30+  =  K4V*D20   '  (153) 

The  ratios  for  K  and  K  are  given  by 

vk2  ■  exp[(-(aG;>Hlo+  -  ag;,dh,o+»/rt!  <im) 

3  2 

where  AG   is  the  standard  free  energy  of  solvation  in  water  of  the  proton  or 
deuteron  in  the  form  of  the  various  oxonium  ions  (18)  . 

For  the  isotopic  species,  differences  of  free  energy  of  solvation  are 

primarily  due  to  the  differences  of  zero  point  energies  of  the  ions,  so  that 

+      + 
differences  of  free  energy  of  solvation  of  H  and  D  can  be  identified  approxi- 
mately with  differences  in  corresponding  heats  of  solvation  calculated  above. 
For  example,  in  0.1N  acid  solution,  containing  H„0,  HDO  and  D90,  the  activities 
of  the  oxonium  ions  are 

a   +  =  10    (mole/liter) 

a   Q+  =  3.8  x  10"4  (mole/liter)  (155) 

—8 
and  a    _+  is  of  the  order  of  magnitude  10   and  a   _+  is  of  the  order  of 

—  1  9 

magnitude  10    (mole/liter)  so  that  only  H~0  and'  DH„0  need  be  considered  as 
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sources  of  protons  or  deuterons  in  acid  solutions. 

The  two  reactions  whose  rates  must  be  compared,  therefore,  are 

Ho0+  +  e~  +  M  ■>   MH  ,  +   H„0  (156) 

3  ads     2 

and  DHo0+  +  e"  +  M  -»•  MD    +  Ho0   .  (157) 

2  ads     2 

Since  thermodynamic  quantities  other  than  the  zero  point  energies  enter  into 

the  calculation  of  the  different  heats  of  activation  for  proton  and  deuteron 

discharge,  it  is  therefore  convenient  to  refer  the  total  heat  content  of 

reactants  and  of  products  to  that  of  the  reference  state,  1/2  H  +  HO  +  M 

in  the  case  of  reaction  (156)  and  to  that  of  1/2  D„  +  H„0  +  M  in  the  case  of 

reaction  (157)  through  the  Born-Haber  cycles. 
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Cycle  I     1/2  H0  +  Ho0  +  M  ■>  H  +  Ho0  +  M  1/2  D° 

2     2  2  n 

H  +  Ho0  +  M  -*■  H+  +  e"  +  Ho0  +  M         I„ 
2  Z  ti 

H+  +  e"  +  Ho0  +  M  +  H„0+  +  e"  +  M       AH  u+  „  ft  -  <j)M 
2         3      M  sHT,H  0    M 

1/2  H0  +  H„0  +  M  -*  Ho0+  +  e"  +  M 
2     2         3      M 

Reference  state  ■*■  Initial  state  of  reaction  (156) 
(AH1>H  "  1/2  DH2  +  XH  +  AHsH+,H20  "  *M 


Cycle  II     1/2  H0  +  Ho0  +  M  -»•  Ho0  +  MH 

2    2         2      ads 

Reference  state  ■>  Final  state  of  reaction  (156) 
(AHf>H  "  1/2  DH2  +  AHadsH 


An  alternative  representation  of  the  above  two  cycles  is 
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[H  +  e  +  HO  +  M] 


Cycle  I 


H 


[H  +  HO  +  M] 


1/2  D° 


[1/2  H2  +  H20  +  M] 


(AHi>H 


AHsH+,H20  "  ^M 


>  [H30  +  eM  +  M] 


Reference  state 


Initial  state  of  (156) 


and 


1/2  D 


[1/2  H2  +  H20  +  M] 


[H  +  HO  +  M] 


(AHf)H 


adsH 


Cycle  II 


->  [Ho0  +  MH    ] 
2      ads 


Reference  state 


Final  state  of  (156) 


For  deuterium  entities  the  Born-Haber  cycles  are 
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[D+  +  e"  +  HO  +  M] 


[D  +  H20  +  M] 


1/2  D 


[1/2  DH  +  HO  +  M] 


(AH.)D 


Cycle  I 


iHsD+,H20  "  *M 


>  [DHo0  +  e„  -I-  M] 
2      M 


Reference  state 


Initial  state  of  reaction  (157) 


1/2  DDH 


[1/2  DH  +  HO  +  M] 


[D  +  HO  +  M] 


(aVd 


Cycle  II 


adsD 


■>  [MD     +  HO] 

ads     2 
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In  these  cycles,  AH  and  AH  are  the  heat  content  changes  in  the  formation  of 

the  initial  and  final  states  from  the  reference  state;  AH  „+  „  „  is  the  heat 

sH  ,H„0 

of  solvation  of  H   in  HO,  AH     is  the  heat  of  adsorption  of  H  at  the  metal 

/     adsn 

surface  and  <$>      is  the  work  function  of  the  metal.   The  values  of  the  terms  are 
given  in  Table  II. 

As  shown  in  Table  II  heat  of  solvation  of  H   in  HO  is  -263  ±  13  Kcal/ 
mole.   The  heat  of  solvation  of  D  can,  therefore,  be  taken  as  -263  -  1.8  (see 
section  above) . 

Using  the  data  given  in  Table  II 

(AH.)R  =  101.6  -  <j.M  ;  (AH.)D  =  100.3  -  ^  .  (158) 

As  seen  (AH.)  values  for  H  and  D  differ  by  1.3  Kcal/mole,  but  this  is 
also  the  difference  between  the  zero  point  energies  of  the  OH  and  OD  bonds 
in  HO  and  H  DO   ions  in  the  initial  state  of  the  reaction.   The  minima  of 
the  potential  energy  curves  for  proton  and  deutron  transfer  from  H„0  or  DH_0 
ions,  respectively,  therefore  lie  at  the  same  potential  energy.   The  minima 
for  the  final  state  curves  differ  in  energy  by  the  quantities  AQ  (see  Appendix 
D) .   The  difference  in  heats  of  activation  for  proton  and  deuteron  transfer 
would,  therefore,  be  equal  to  the  difference  in  AH.  values  for  H  and  D,  plus 

£AQ,  where  3  =  0.5  is  the  symmetry  factor.   An  additional  effect  arises  due  to 

+      + 
the  slightly  different  shapes  of  the  Morse  curves  for  H  and  D  dissociation 

from  HO  and  H  DO  ,  respectively.   Evaluation,  as  described  by  Parsons  and 

Bockris,  (28),  of  the  Morse  constants  "A"  for  these  dissociation  processes, 

gives 

AH+H   =  1.62  x  108  and  A^  Q   =  1.65  x  108cm-1  .  (159) 


Substituting  the  above  values  of  "A"  into  the  Morse  equation  using  the  heats 
of  solvation  and  zero  point  energies  given  in  Table  II  for  the  H  -  H„0  and 


Table  II.   Energy  quantities  used  in  the  calculations 
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Quantity 


Value  in  Kcal/mole 


Reference 


Ionization  potential  of  H 
or  D 

Bond  dissociation  energy 
D°  of 

H2 
HD 

D2 

Proton  affinity  of  water 
Deuteron  affinity  of  water 


Total  heat  of  solvation  of 

i2c 


H+  in  H„0 


Total  heat  of  solvation  of 
D+  in  H„0 


Total  heat  of  solvation  of 
D+  in  DO 

OH  bond  dissociation  energy 
in  HO 

OD  bond  dissociation  energy 
in  HOD 


Heat  of  ionization  of  H„0 
in  HO 

Heat  of  ionization  of  HOD 
in  H20  (to  give  H  D0+  +  0H~) 

Heat  of  ionization  of  D„0 
in  D90 

Bond  dissociation  energy  of 

CuH,  CuD 


I„  =  313.00;  I  =  313.09 
n  D 


Du0  =  103.12 

H 


DHD°  =  103.94 


D   =  104.92 
2 


-182 
-183.8 
-263  *  13 

-264.8 

-264.2 

120 

121.5 

13.50 

13.89 

14.50 


(21) 


(26) 
(26) 
(26) 


(27,28) 
(29) 

(27,28) 

(27,28) 

(30) 

(30) 

(31) 

(18) 

(27) 


DCuH=66-8'  DCuD=67-6 


(32) 


Table  II  (continued) 
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Quantity 


Value  in  Kcal/mole 


Reference 


AgH,  AgD 
PbH,  PbD 
HgH,  HgD 
N  in- 
Zero  point  energy  of 
CuH,  CuD 

AgH,  AgD 

PbH,  PbD 

HgH,  HgD 

NiH,  NiD 

Zero  point  energy  of  OH 
bond  in  H  0+ 

Zero  point  energy  of  0D 
bond  in  H  D0+ 


D°   =  57.7;  D°   =  58.4 
AgH  AgD 


D;bH  -  37-°-  D;bD  -  37-6 


D°    =  57.0-  D°    =  57  5 
UHgH    3/,U'   HgD    3/-D 


DNiH  =  71'° 


ECuH  "  2'78'  ECuD  =  1'99 


E°   =  2.53;  E°   =  1.79 
AgH  AgD 


E;bH  =  2-20'  E;bD  =  i-7° 


E°   =  1.99;  E°   =  1.43 
Hgh  HgD 


E°  „  =  2.27;  E°   =  1.61 
NiH  NiD 


5.0 


3.73 


(33) 
(33) 
(28) 
(33) 


(33) 
(28) 

(28) 
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+  + 

D   -  HO  interaction,  show  that  the  Morse  function  for  the  dissociation  of  D 

from  HO  rises  more  steeply  than  that  for  H  .   Potential  energy  diagrams  show 

that  in  the  activated  state  of  the  discharge  reaction,  the  proton  or  deuteron 

bond  to  oxygen  has  to  be  stretched  by  about  0.2A°  from  the  equilibrium  separa- 

+  + 

tion.   At  this  distance,  the  Morse  curve  for  the  dissociation  of  D   from  DH„0 

is  0.6  Kcal  higher  than  that  for  H+  from  H  0+  (18).  With  3  =  0.5  the  activa- 
tion energy  for  D  transfer  will  hence  be  about  0.3  Kcal  higher  than  that  for 
H  on  this  account.   The  total  difference  of  heat  of  activation  is  therefore 

(AH^  -  AH„)  =  1.3  +  1/2(0.36)  +  1/2(0.60)  =  1.78  Kcal/mole.       (160) 
D      H 

Using  the  ratio  C  /C__  =  1/5500  in  ordinary  water  solutions;  a   +/a    +  = 

D   H  H„U    DH~U 

—  1  — / 

10   /3.8  x  10   as  shown  above  in  Eq.  (155)  and  the  fact  that  the  probability 
of  proton  being  discharged  from  HO   is  3  times  that  of  a  deuteron  from  DH  0 
in  the  double  layer  (34),  the  value  of  S  is  given  by 

S  -  (1/5500) (10_1/3. 8  •  104)  3  exp(1.78/RT) 

=  2.9  .       (161) 

(b)   Discharge  in  water  containing  either  pure  HO  and  protons  or  pure  D_0 
and  deuterons. 

Here  the  relative  ratios  of  the  simple  discharge  reaction,  from  H„0   in 
HO  and  from  DO   in  DO  for  the  same  ionic  concentrations,  are  compared. 

H„0+  +  e  +  M  ->  Ho0  +  MH  ,  I 

3  2      ads 

Do0+  +  e  +  M  ->  Do0  +  MD  II 

3  2      ads 

1/2  H„  +  H^O  +  M  is  taken  as  the  reference  state  for  reaction  I  and  1/2  D„  + 

Do0  +  M  as  the  reference  rate  for  reaction  II.   For  reaction  I  (AH.)   = 

i.  In 

101.56  -  d>   as  calculated  above.   For  reaction  II,  Born-Haber  cycles  are 

M  J 


D  +  DO  +  M 


42 


1/2  D 


1/2  D  +  DO  +  M 


iHsD+,D20  "  *M 


(AH.)D  -  1/2(DJ)  +  ID  +  AH^     -  ^ 


=  1/2(104.92)  +  313.09  -  2.64.2  - 


M 


=  101.35  -  $ 


M 


The  difference  (101.56  -  101.35)  =  0.21  is  not  equal  to  the  difference  in  zero 
point  energy.   Therefore  the  minima  of  the  initial  state  potential  energy 
curves  differ  by  the  difference  of  zero  point  of  OH  bond  in  H„0  and  OD  bond 
in  DO  .   The  minima  of  the  final  state  potential  energy  curves  differ  by  0.36 
(Table  III).   Hence 

(AH*  -  AH*)  =  0.21  +  1/2[(1.27  -  .21)  +  (0.56  -  0.2)]  +  1/2(0.60) 
D  n 


=  1.22  Kcals/mole 


(162) 


The  activity  ratio  factor  a  n+/a   +  and  the  partition  function  ratios  are 
both  unity.   Hence  S  from  Eq.  (135)  is  given  by 
S  =  exp((1.2  x  1000)/(1.987  x  298)) 

=  7.7  (163) 
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Table  III.   Zero  point  energies  and  (AHf)  values  for  the 
simple  discharge  reaction  at  various  metals 


Difference 

AQ  of  depth 

Zero  point  energy  of  MH/MD  p.e. 

E  of  MH  or  MD         AE_  ,„,  ._         curve  minima 

o 0  MH-MD 

1.99  0.56  0.36 

1.43 

2.78  0.79  0.39 

1.99 

2.53  0.74  0.44 

1.79 

2.20  0.64  0.44 

1.60 


[H/MD 

(AHf) 

HgH 

-5.4 

HgD 

-5.6 

CuH 

-15.2 

CuD 

-15.6 

AgH 

-6.1 

AgD 

-6.4 

PbH 

14.6 

PbD 

14.4 
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Case  II:   Electrochemical  Desorption 

(a)   Discharge  in  water  containing  the  normal  isotopic  mixture  of  H  and 
D. 

The  following  reactions  are  possible. 


H„0  +  e  +  MH    +  Hn  +  M  +  H„0 
3  ads     2        2 


DH_0+  +  e  +  MH  J   -*  HD  +  M  +  Ho0 
2  ads  2 


II 


H„0+  +e  +  MD    -*  HD  +  M  +  H^O 
3  ads  2 


III 


Here  the  heat  contents  of  states  in  the  reactions  are  referred  to  that  of 
H  +  HO  +  M  or  HD  +  HO  +  M.   The  Born-Haber  cycle  for  the  initial  state  of 
reaction  I  is 


H    +  e  +  H    +M+  Ho0 
gas        gas        2 


11 


2H    +  M  +  Ho0 
gas 


H0     +  M  +  H_0 
2  gas        2 


(AH.) 


AH  u+  u  n  +  AH  a     u  "  ^u 

sHT,H  0     adsH    MH 


>  Ho0  +  MH    +  «» 
3       ads    M 


(AVh  =  DH2  +  h   +  aHsH+,H20  +  AHadsH  "  *MH 


Substituting  the  values  for  example  for  Ag 
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=  103.12  +  313  +  (-263)  +  (-57.7)  -  <|> 

=  95.42  -  $-_-   .  (164) 

MH 

For  a  similar  cycle  for  reaction  II 

(AVd  via  DH20+  "  °HD  +  h  +   AHsD+,H20  +  AHadsH  "  *MH 

=  94.53  -  $._,   .  (165) 

MH 

and  analogously  for  reaction  III 

(AVD  via  MD  „   "  DHD  +  h   +  AHsH+,H  0  +  AHadsD  "  *MH 

ads  2 

=  95.54  -  <j>__.   .  (166) 

MH 

The  heat  contents  of  the  final  states  are  those  of  physically  adsorbed  H„  at 
the  metal  for  reaction  I,  or  HD  for  II  and  III,  and  do  not  differ  significantly, 
However,  the  minima  of  the  final  state  potential  energy  curves  for  the  forma- 
tion of  H  and  HD  will  differ  by  the  difference  of  the  zero  point  energies  of 
H  and  HD,  i.e.,  by  0.82  Kcal/mole  which  is  the  difference  in  dissociation 
energies  of  H?  and  HD.   With  6  =  0.5,  this  gives  a  contribution  of  0.41  Kcal/ 
mole  to  the  difference  of  heat  of  activation  for  H  and  HD  liberation  by  the 
electrochemical  desorption  mechanism.   Similarly,  for  reactions  I  and  II,  the 
difference  of  the  minima  of  the  initial  state  curves  will  be  equal  to  the 

difference  of  heat  contents  for  the  two  initial  states  (namely,  95.42  -  94.53 

+       + 
Kcal/mole)  minus  the  difference  of  zero  point  energies  of  the  OH  and  OD  bonds 

(namely,  1.27  Kcal/mole)  i.e.,  0.38  Kcal/mole,  the  curve  for  the  D   transfer 

being  the  higher.   Including  the  effect  (0.6  3Kcal)  due  to  the  difference  of 

Morse  functions  for  H  and  D  dissociation,  the  total  difference  of  heat  of 

activation  between  reactions  I  and  II  will  be 

95.42  -  94.53  +  1/2(0.82)  +  1/2(0.38)  +  1/2(0.6)  =  1.8  Kcal/mole   . 
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With  aH  0+/aH  D0+»  the  ratio,  r,  of  rates  of  H   and  HD  production  by  reactions 
I  and  II  respectively,  is 

r  =  2.64  x  102  exp(1800/591)  -  1.68  x  104  -  v   /v  (167) 

where  v   and  v   are  the  velocities  of  production  of  H_  and  HD.   The  value  of 
r  is  related  to  separation  factor,  S,  by  relation 

S  "  (CD/CH>soln    (l/2)vHn *  <  VCH>sol»<2r  +  «  ■  <168) 

tiu 

Since  in  the  production  of  HD,  one  H  atom  accompanies  each  D  atom  liberated. 
If  v    >>v   ,  as  in  the  case  considered  above 

2 

S  =  (C_/C„)   . =  (C_/Cu)   .   2r 

D  H  soln  v       D   H  soln 
HD 

=  (1/5500)  x  2  x  1.68  x  104  =  6.1  (169) 

i.e.  almost  twice  the  value  arising  if  simple  discharge  were  rate  determining. 

In  the  calculation  of  S  it  has  been  assumed  that  the  ratio  of  H  .   to  D     is 

ads     ads 

equal  to  the  ratio  of  activities  of  H„0   and  H  DO   in  the  solution.   When  the 

electrochemical  desorption  (III)  is  rate  determining  but  faster  than 

MH    +  MH    ■*  2M  +  H0  , 
ads     ads         2 

the  steady  state  concentration  of  adsorbed  H  has  been  calculated  by  Conway  (18) 

and  the  value  of  S  then  calculated  for  this  mechanism  was  0.7,  i.e.,  a  value 

not  observed  experimentally.   Hence  atomic  desorption  is  an  unlikely  path  in  HD 

evolution.   However,  the  conditions  for  which  S  would  be  0.7,  are  limiting  and 

correspond  to  those  for  which  Tafel  slope  of  0.038  for  the  electrochemical 

desorption  mechanism  is  predicted  (35),  i.e.,  almost  twice  that  arising  if 

simple  discharge  were  rate  determining. 

(b)   Discharge  of  HO   and  DO   ions  in  pure  H„0  and  D»0  respectively  at 
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the  same  acid  concentration. 

For  this  case,  the  reactions  are 


Ho0+  +  e  +  MH  j   ■*  h\  +  h\0  +  M 
3  ads     2     2 


Do0  +  e  +  MD  ,  ->   D0  +  Do0  +  M   .  IV 

3  ads     2     2 


(AH.)   for  reaction  I,  at  silver,  is  95.42  -  <J>    Kcal/mole  and  for  reaction 
1  ri  A.gH 

IV,  (AH.K  =  104.92  +  313.09  -  58.4  -  264.2  =  95.41  -  <f»A  ^.   The  zero  point 
l  D  AgD 

energy  levels  for  the  two  initial  states  are  therefore  coincident,  so  that  the 
minima  of  the  initial  state  potential  energy  curves  will  differ  by  1.27  Kcal/ 
mole,  that  for  DO   being  higher.   The  minima  of  the  final  state  potential 
energy  curves  will  differ  by  the  difference  of  zero  point  energy  of  H„  and  D~ , 
namely,  1.80  Kcal/mole.   The  heats  of  activation  hence  differ  by 

AH*  -  AH*  =  (95.42  -  95.41)  +  1/2(1.27)  +  1/2(1.80) 

D  n 

=  1.545  Kcal/mole   .  (170) 

The  activity  ratio  a   +/a        +   is  unity  and  the  partition  function  ratio  in 

H3°   D3° 
Eq.  (135)  is  2  for  the  above  case  (18). 

The  hypothetical  separation  factor  arising  under  the  limiting  conditions 

of  acid  solutions  in  pure  H90  and  D?0  is,  therefore 

S  =  2  exp[1545/RT]  =  27  .  (171) 

Case  III:   Atomic  Desorption  (or  Atomic  Recombination) 

(a)   Discharge  in  water  containing  the  normal  isotopic  species  of  H  and  D, 

The  reactions  in  ordinary  water  solutions  in  this  case  are 

MH    +  MH  J   -+  2M  +  H„  I 

ads     ads         2 

MH    +  MD    -»■  2M  +  HD  II 

ads     ads 

At  high  surface  coverage  the  difference  of  activation  energy  for  reactions  I 
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and  II  is  determined  by  the  difference  between  the  zero  point  energies  of  one 

MH    and  one  MD    bond.   No  data  for  the  vibrational  frequency  of  PtH  (for 
ads  ads  n    ' 

which  it  is  known  that  the  above  reaction  is  rate  determining  (18))  is  avail- 

13 
able,  but  it  is  probably  close  to  that  for  NiH,  namely  v  „.,,  =  4.75  x  10 

J      0N1H 

and  v  M._  =  3.36  x  10   /sec. 
0N1D 

The  partition  function  ratio  is  4/3  (18)  so  that  the  ratio  of  the  rate 

of  evolution  of  H   to  that  of  HD  is 

r  -  (4/3)exp[l/2((hvo)N.H  -  1/2  (hv^^) /kT]  =  4.1         (172) 

-27 
where  h  =  Planck's  constant  =  6.625  x  10    ergs 

k  =  Boltzmann  constant  =  1.38  x  10    ergs/deg 
and  since  S  =  (2r  +  1)  as  shown  in  Eq .  (168),  therefore  S  =  9.2.   On  an  acti- 
vated platinum  surface  in  1.5  N  sulphuric  acid,  the  experimental  value  of  S 
is  6.4  (36,  37).   Lower  values  have  also  been  reported  and  it  has  been 
suggested  that  these  arise  not  from  an  electrolytic  separation,  but  because 
the  HDO  +  H„  J  HD  +  H_0  equilibrium  is  set  up  at  the  catalytically  active 
metal;  this  is  disproved,  however,  by  the  work  of  Okamoto  (38).   The  value 
calculated  is  an  upper  limit,  since  the  heat  of  adsorption  of  H  is  diminished 
as  the  coverage  of  the  metal  by  H  approaches  saturation,  so  that  the  zero 
point  energies  of  the  MH  bonds  would  be  expected  to  be  lowered  proportionally. 
Smaller  values  of  S  would  hence  arise;  the  minimum  value,  when  the  exponential 
in  Eq.  (172)  is  unity,  is  2.7. 

(b)   When  atomic  desorption  reaction  takes  place  in  pure  H_0  and  DO  the 
reactions  to  be  considered  are 

MH  J   +  MH     +  M  +  H0  I 

ads     ads        2 

MD    +  MD    ■*■   M  +  D0  III 

ads     ads        2 

The  ratio  of  partition  functions  is  2  (18) .   Hence 
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S-  2  exp[(l/2(hvo)NiH-  l/2(hvo)N.D)] 

=  6   .  (173) 

Results  of  the  derivations  presented  above  for  various  reaction  mechanisms 
are  summarized  in  Table  IV. 
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2.5  Determination  of  Corrosion  Rates  from  Overpotential  Measurements 

Stern  and  Geary  (34)  showed  that  for  corroding  systems  controlled  by 
activation  polarization  a  linear  relationship  is  expected  in  the  region  where 
the  polarized  potential  is  close  to  the  corrosion  potential.   For  these  condi- 
tions, the  following  equation  was  derived. 

b  b 

lAiJ       (2.3)(i     )(b   +  b  )  U  } 

E+o  corr   c    a 

where  AE/Ai  is  the  polarization  resistance,  the  constants  b   and  b   are  the 
'  r  '  c      a 

slopes  of  the  logarithmic  local  cathodic  and  anodic  polarization  curves,  and 

i     is  the  corrosion  current, 
corr 

Mueller  (40)  derived  the  polarization  curve  of  alloys  from  those  of  the 
single  metallic  components.   According  to  this  author  the  rate  of  dissolution 
of  an  alloy  in  the  active  state  approaches  that  of  the  main  metallic  component 
of  lowest  rate  of  dissolution  within  certain  limits. 

Usually,  the  conditions  in  an  alloy  are  complicated  as  the  alloy  might 
not  represent  a  mixture  of  the  crystals  of  pure  metals.   Each  metal  phase 
might  contain  a  certain  percentage  of  all  the  metals  present  or  all  metals 
might  be  contained  in  a  mixed  crystal  as  in  the  case  of  stainless  steel.   If 
the  free  energy  of  the  metal  components  is  not  essentially  changed  by  alloying, 
derivation  of  the  polarization  curve  of  the  alloy  from  those  of  the  metallic 
components  might  be  possible.   The  condition  of  the  absence  of  an  essential 
drop  of  the  free  energy  by  alloying  is  the  opposite  of  that  of  the  validity 
of  Uhlig's  electron  configuration  theory  (41),  according  to  which  the  d-electron 
vacancies  in  one  of  the  component  metals  are  filled  by  electrons  originating 
from  some  other  metallic  component  or  hydrogen.   This  mechanism  would  involve 
a  drop  of  free  energy  (42) . 
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The  conditions  of  the  applicability  of  Mueller's  theory  to  stainless 
steels  are  satisfied,  because,  according  to  Evans  (42),  the  free  energy  of  Fe, 
Cr,  and  Ni  is  not  essentially  changed  by  alloying  to  form  stainless  steels. 

On  the  basis  of  the  above  theory,  the  overall  current  density  of  the 
alloy  may  be  written  as 


AAA 
x    .    y    .    ; 

a    x  A     y  A     z  A 


i  =i  T^+i  T^+i  IT  (175) 


where  i  =  current  density  of  the  alloy 

a  J  J 

i  ,  i  ,  i  =  current  densities  of  the  single  metals  X,  Y,  Z 
x   y    z 

A  ,  A  ,  A  =  exposed  areas  of  the  metals  X,  Y,  Z 
x   y    z 

A  =  total  exposed  area  of  the  alloy 

If  the  concentrations  of  the  single  metals  on  the  surface  are  not 
different  from  those  in  the  bulk,  the  ratio  of  exposed  areas  of  the  single 
metals  to  total  exposed  area  of  the  alloy  would  equal  the  ratio  of  corre- 
sponding concentrations  of  the  metals  in  the  alloy  to  the  specific  gravity  of 
the  respective  metals.   Hence  Eq .  (175)  can  be  written  as 

C       C       C 

i  -  i  -*  +  i  -^  +   i  ^  (176) 

a    x  S     y  S     a  S 
x       y       z 

where  C  ,  C  ,  C  =  concentrations  of  single  metals  x,  y  and  z  expressed  in 
x   y   z 

weight  per  unit'  volume  of  the  alloy 
S  ,  S  ,  S  =  specific  gravity  of  the  single  metals  X,  Y,  Z. 

The  rate  at  which  metal  corrodes  is  proportional  to  i    ,  the  corrosion 

r   r  corr 

current,  defined  by  (5) 

i 

E     =  -b  log  COrr  (177) 

corr  l  \    >  •  j 


where  E     is  the  corrosion  potential  (see  fie.  below) 
corr  r  ° 
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Eq.  (177)  may  be  used  to  calculate  the  corrosion  current  if  b,  i  and  the 

corrosion  potential,  E    ,  are  known.   The  slope  b  is  determined  graphically 

corr 

from  the  Tafel  line  of  the  cathodic  polarization  curve.   The  exchange  current 
is  determined  by  extrapolation  of  the  Tafel  region  to  the  reverse  hydrogen 
potential  and  the  corrosion  potential  from  the  cathodic  and  anodic  polariza- 
tion curves  in  terms  of  i     for  alloys  as 

corr 


C  C  C 

X  V  z 

corr  a     corr  x  S      corr  y  S      corr  z  S 

x  y  z 


(178) 


Eq.  (178)  indicates  that  the  rate  at  which  an  alloy  corrodes  is  given  by 
summing  the  corrosion  rates  of  individual  elements  multiplied  by  its  fraction 
of  surface  area  exposed.   This  is  true  for  a  freshly  polished  surface.   After 
the  corrosion  has  begun,  the  situation  is  much  more  complex  than  indicated 
in  Eq.  (178),  especially  in  case  of  passive  alloys. 
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3.0   EXPERIMENTAL 

3.1   Description  of  Equipment 

Fig.  4  shows  the  equipment  set-up  used  to  obtain  the  overvoltage  data 
for  HC1  in  H20  and  DC1  in  D20  on  stainless  steel  and  nickel  electrodes.   The 
system  consisted  of 

A.  Electrolytic  Cell 

B.  Hydrogen/Deuterium  Purification  Unit 

C.  Water/Heavy  Water  Purification  Unit 

A.   Electrolytic  Cell 

The  electrolytic  cell,  shown  in  Fig.  5  was  made  of  Pyrex  glass  and  was 
comprised  of  the  following  three  compartments. 

Reference  Electrode,  I 

Cathode,  II 

Anode,  III 

I   Reference  Electrode  Compartment  (4  cm  diam,  12  cm  high)  consisted  of 

2,  hydrogen  electrode,  with  1  cm  x  1-1/2  cm  platinized  platinum 
electrode,  1. 

3,  $  45/50  ground  joint  cover. 

4,  water  trap,  6  cm  long  by  2  cm  diam. 

II   Cathode  Compartment  (6-1/2  cm  diam,  22  cm  high) 

5,  test  electrodes. 

6,  thermo-well  made  of  glass  for  holding  a  thermometer. 

7,  $  29/42  ground  glass  joint  with  water  seal. 

8,  $  71/60  ground  glass  joint  with  water  seal. 

9,  four  syringe  type  ground  glass  joints  to  provide  a  leakproof 
adjustable  arrangement. 
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10,  water  trap,  6  cm  by  2  cm. 

Ill  Anode  Compartment  (6-1/2  cm  diam,  22  cm  high) 

11,  fritted  glass  disk,  1  cm  diam,  fitted  at  the  bottom  of  1-1/2  cm 
diam  tubing  which  carries  the  anode. 

12,  platinized  platinum  electrode,  1  cm  by  1-1/2  cm. 

13,  two  1  cm  by  1  cm  platinized  platinum  electrodes  for  measurement 
of  conductivity. 

14,  $  19/38  ground  glass  joint  with  water  seal,  carries  the  anode,  12. 
15  and  31,  water  trap,  6  cm  by  2  cm. 

16,  $  71/60  ground  glass  joint  with  water  seal. 

17,  Luggin  capillary:  a  tapered  tubing  1-1/2  cm  long  protruding 
into  the  cathode  compartment  and  having  1  mm  capillary  mouth. 

18,  $  10/30  ground  glass  joint  with  water  seal  for  transmitting  low 
conductivity  water. 

19,  $  12/30  ground  glass  joint  on  which  a  trident  shaped  graphite 
electrode  could  be  fitted  for  pre-electrolysis . 

20  and  26,  $  10/30  ground  glass  joints  for  hydrogen,  deuterium  or 
nitrogen  supply  to  the  anode. 

21,  $  10/30  ground  glass  joint  used  as  inlet  to  the  anode  compartment. 

22,  drain. 

25,  main  gas  inlet  $  10/30  ground  glass  joint,  water  seal  type. 

23,  24,  27,  28,  29,  and  30,  glass  valves. 

B.   Hydrogen/Deuterium  Purification  Unit,  Fig.  6 

Commercial  hydrogen  gas  usually  contains  trace  impurities  of  oxygen, 
carbon  monoxide,  carbon  dioxide,  hydrocarbons  and  water  vapor.   To  get  rid  of 
these  impurities,  H»  was  first  passed  through  a  series  of  glass  bubblers,  each 
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approximately  3  cm  diam  by  16  cm  long;  the  first  two  containing  dry  silica  gel 
to  remove  moisture;  the  next  two  containing  hopcalite  to  remove  carbon  monoxide 
in  absence  of  moisture.   The  hydrogen  gas  was  then  passed  through  palladized 
asbestos  (5%  Palladium)  tube  immersed  in  an  electrically  heated  oven  at  500°C 
to  remove  oxygen  by  catalytic  combination  with  hydrogen.   Further, bubblers  of 
silica  gel  were  used  to  remove  water  produced  in  the  palladium  asbestos  tube. 
Four  glass  traps  followed;  two  containing  silica  gel  and  activated  charcoal, 
to  remove  moisture  and  some  gaseous  impurities,  and  two  were  kept  empty.   The 
last  two  traps  were  immersed  in  dry  ice  to  trap  any  trace  of  moisture  escaping 
adsorption  by  the  silica  gel. 

The  same  arrangement  was  used  for  purification  of  deuterium  gas.   All  the 
tubes  containing  silica  gel,  hopcalite,  palladized  asbestos  and  activated 
charcoal  were  dried  in  the  oven  at  250°C  for  about  48  hours  to  ensure  complete 
removal  of  moisture. 

C.   Water  Purification  Unit 

Water  used  to  prepare  solutions  for  electrochemical  measurements  was 
purified  by  redistillation  of  predistilled  water  in  an  all  Pyrex  glass  appa- 
ratus.  About  60%  of  the  distillate  was  used.   Its  specific  conductance  was 
1  x  10   mhos/cm.   Purified  hydrogen  gas  was  passed  all  the  time  during 
distillation.   H„  gas  passed  was  purified  according  to  the  method  described 
above. 

Electrical  Set-Up  (Fig.  7) 

The  electric  circuit  used  consisted  of: 

1,  potentiometer:  Type  K-3  Universal,  Guarded,  made  by  Leeds  &  Northrup. 

2,  storage  cell:  200  ampere  hours;  2  volts. 

3,  standard  cell:  Epply  unsaturated*  internal  resistance  less  than 
500  ohms. 
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4,  galvanometer:  sensitive  reflecting  type,  used  as  a  null  point 
indicator, 

5,  electrometer:  Keithley  Model  600A  Electrometer,  made  by  Keithley 

-13 
Instruments;  a  direct  reading  ammeter  from  3  amperes  to  10    amperes 

full  scale.   The  current  could  be  measured  with  an  accuracy  of  3%  of 

-13 
full  scale  from  3  to  10    amperes. 

6,  condenser:  2.5  microfarads,  connected  in  parallel  with  the  input  leads 
of  the  electrometer.   This  arrangement  avoided  the  effect  of  extran- 
eous currents  on  the  measurements. 

7,  pH-meter:  Beckman  Model  G,  used  to  serve  the  dual  purpose  of:  (i) 
measuring  the  pH  of  the  solution  and  (ii)  acting  as  a  null  point 
indicator  for  the  measurement  of  potential  differences. 

8,  power  supply:  Heathkit  variable  voltage  type;  used  for  pre-electrolysis 
and  platinization  of  electrodes. 

9,  variable  resistance  box;  15  ohms  to  10  megohms. 

10,  6  dry  cells,  l.J  volts  each. 

11,  conductivity  bridge:  Model  RC,  made  by  Industrial  Instruments. 

12,  thermostat:  for  studies  on  temperature  dependence;  electronically 
controlled  to  keep  constant  temperatures  within  ±   0.1°C. 

3.2  Techniques 

Cathode  Preparations 

Stainless  steel  electrodes  used  as  cathodes  were  of  the  stainless  steel- 
304  type,  1/16"  diam  wire,  supplied  by  United  Steel  Corporation,  of  known 
composition  (see  Appendix  B) .   The  sample  was  already  annealed  by  the  suppliers 
and  hence  further  annealing  prior  to  electrochemical  measurements  became  un- 
necessary.  Each  stainless  steel  cathode  consisted  of  1-1/4"  of  the  above 
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specimen  attached  to  a  short  length  of  28  SWG  Pt  by  spot  welding.   The  Pt 
wire  was  then  fused  into  the  Pyrex  glass  tubing  carrying  the  electrode. 

Each  electrode  was  first  polished  with  an  emery  paper,  then  electro- 
polished  in  a  solution  consisting  of  65%  by  volume  H_SO  ;  22.5%  H?0  and 

12.5%  H„PO, .   The  electrolyte  temperature  was  maintained  at  80°C  and  a 

2 
current  of  50  ma/cm  was  passed  for  ten  minutes.   The  electrode  was  then 

washed  thoroughly  with  distilled  water  and  immediately  fixed  in  the  cathode 

compartment  where  pure  hydrogen  gas  atmosphere  was  maintained. 

Nickel  cathodes,  on  the  other  hand,  were  made  from  a  spectroscopically 

pure  sample  of  known  analysis  supplied  by  Johnson,  Matthey  &  Co . ,  Ltd. (see 

Appendix  C) .   Each  piece  (1  cm  x  2  cm,  1/16"  thick)  was  spot  welded  to  Pt  and 

sealed  into  glass.   The  Ni  cathodes  were  uniformly  polished  with  emery  paper, 

washed  first  with  acetone  and  then  with  conductance  water  before  being  fixed 

in  the  cathode  compartment. 

Preparation  for  Overvoltage  Measurements  in  HC1-H-0  Mixtures 

The  glass  apparatus  was  cleaned  in  acid  solution,  followed  by  several 
washings  with  water  of  increasing  purity.   The  final  washing  was  done  with 
conductance  water  and  was  carried  out  until  the  conductivity  of  the  washings 
was  less  than  5  x  10   mhos/cm.   The  whole  system  was  then  assembled,  all  the 
traps,  joints  and  valves  were  sealed  with  conductance  water  and  the  whole 
system  was  tested  for  leakage.   The  hydrogen  purification  unit  was  evacuated 
and  flushed  with  nitrogen  several  times  to  ensure  removal  of  air,  in  order  to 
avoid  explosions.   The  temperature  of  the  electrical  oven  surrounding  the 
palladized  asbestos  tube  was  maintained  between  450°C  and  500°C.   Air  was 
expelled  from  the  cell  by  passing  a  fast  current  of  pure  H„. 

A  known  volume  of  HC1  solution  of  predetermined  concentration,  prepared 
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by  mixing  the  appropriate  amount  of  triply  distilled  azeotropic  mixture  of 
HC1-H  0  with  conductance  water,  was  then  transferred  into  the  anode  compart- 
ment of  the  cell.   Transfer  of  conductance  water  from  the  distillation 
apparatus  to  the  mixing  flask  and  of  the  prepared  mixture  to  the  cell  were 
done  out  of  contact  of  air  by  applying  a  pressure  of  purified  H„  gas.   The 
freshly  polished  electrode  was  then  transferred  into  the  cathode  compartment. 
The  concentration  of  the  HC1  solution  was  checked  by  measuring  the  conductiv- 
ity.  The  solution  was  allowed  then  to  be  shared  between  the  anode  and  cathode 
compartments,  the  electrode  being  kept  above  the  level  of  the  solution, 

Pre-electrolysing  the  solution  before  overvoltage  runs  was  found  to  be 
of  extreme  importance  in  obtaining  reproducible  results,  as  shown  below  and 
as  also  indicated  by  Bockris  (7).   For  pre-electrolysis  the  electrical  circuit 

was  connected  as  shown  in  Fig.  7.   The  Pt  electrode  was  dipped  into  the  solu- 

.   2 
tion  and  the  current  adjusted  to  10  ma/cm  .   A  small  current  of  H0  gas  was 

allowed  to  pass  into  the  anode  tube  to  drive  out  the  gasses  produced  during 
the  period  of  pre-electrolysis,   A  period  of  15  hours  of  pre-electrolysis  was 
found  sufficient  to  give  reproducible  results  within  the  experimental  limits. 
Longer  periods  of  pre-electrolysis  did  not  improve  the  reproducibility. 

At  the  end  of  the  pre-electrolysis  period,  the  system  was  ready  for  the 
overvoltage  measurement.   Hydrogen  gas  was  then  allowed  to  pass  in  the  refer- 
ence electrode  compartment  for  15  minutes  before  measurements  were  started. 

-3 
Overpotential  measurements  were  taken  for  current  densities  of  zero  to  10 

amps /cm  . 

Temperature  Dependence  Measurements 

In  temperature  dependence  work  the  thermostat  was  adjusted  to  the  appro- 
priate temperature  a  few  hours  before  the  actual  measurements  were  taken  to 
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make  sure  that  the  inside  of  the  electrolytic  cell  attained  the  desired  tem- 
perature.  Temperatures  below  ambient  were  obtained  by  fixing  an  ice  box 
inside  the  thermostat.   Measurements  were  taken  in  this  work  at  15°,  25°, 
35°,  45°  and  55°C. 

Measurements  in  DC1-D-0  Mixtures 

The  general  experimental  procedure  for  conducting  runs  in  D_0  solutions 
was  with  a  few  exceptions,  very  similar  to  that  used  for  ordinary  water 
solutions.   An  identical  but  smaller  electrolytic  cell  was  used  to  conserve 
as  much  of  D„0  as  possible.   Deuterium  chloride  solutions  were  prepared  by 
bubbling  pure  DC1  gas  supplied  by  Bio-Rad  Laboratories  through  known  volume 
of  D_0.   A  small  amount  of  the  solution  was  titrated  against  standard  sodium 
carbonate  solution,  and  a  calculated  volume  of  D„0  was  added  to  adjust  the 
solution  to  the  required  normality. 

D„  gas  was  purified  in  the  same  manner  as  H„ .   Extreme  care  was  taken  to 
make  the  system  leakproof .   At  the  beginning  of  each  run  the  whole  system  was 
completely  flushed  with  D_  gas  and  during  the  course  of  the  run  D~  gas  was 
circulated  between  the  electrolytic  cell  and  the  purification  system  by  a 
Varistaltic  pump  to  save  on  the  amount  consumed. 

Special  Observations 

It  should  be  emphasized  here  that  unless  extreme  care  is  taken  before 
and  during  the  runs  it  would  be  very  difficult  to  obtain  reproducible  results 
The  following  points,  however,  were  very  helpful  in  aquiring  reproducibility. 
1.   Long  pre-electrolysis  to  purify  the  solution  from  trace  impurities. 
After  trying  different  periods  of  pre-electrolysis  it  was  concluded 
that  the  optimum  period  was  15  hours  for  HC1-H-0  and  26  hours  for 
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DC1-D  0  solutions.   Longer  periods  did  not  improve  the  results  to 

any  noticeable  extent. 
2.   Electropolishing  of  the  electrode  surface:   Some  disagreement 

between  results  obtained  by  different  investigations  in  this  field 

resulted  from  the  irreproducible  finishing  of  electrode  surfaces. 

In  this  work  electropolishing  in  the  solution  described  above 

proved  to  be  very  effective  in  that  direction. 
Although  these  two  points  were  of  very  noticeable  help  yet  it  should  be 
mentioned  that  the  reproducibility  of  runs  under  very  identical  conditions 
was  still  within  limits.   The  statistical  approach  discussed  later  in  Section 
4.1  was,  therefore, acquired  to  take  care  of  this  problem. 
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A.O   RESULTS  AND  DISCUSSION 
4.1  Measurements  on  Ni  Electrodes 

As  indicated  above  in  Section  3.2  it  is  difficult  to  obtain  overpotential 
results  unless  extreme  care  is  taken  in  the  preparation  of  solution  and  elec- 
trode surface.   To  check  the  results  obtained  in  this  work  with  those  obtained 
by  other  investigators,  work  was  carried  on  Ni  electrodes  in  0.01  N  HC1-H„0 
and  0.01  N  DC1-D  0  solutions  which  have  been  previously  studied  by  Bockris 
(43)  and  Conway  (44).   Tables  V  and  VI  show  typical  runs  on  Ni  in  0.01  N  HC1- 
H~0  and  0.01  N  DC1-D?0  solutions.   Because  of  the  difficulties  mentioned  above 
it  became  necessary  to  repeat  measurements  under  identical  conditions  and  treat 
the  results  statistically  as  follows. 

Tafel  slope  b  and  exchange  current  density  i  were  obtained  for  each  n  - 

i   run  by  least  square  fitting  of  the  data  in  the  linear  region.   An  IBM  1410 

program  was  written  for  this  purpose  (see  Appendix  A) .   Variance  of  the  values 

of  i  and  b  obtained  from  least  square  fit  were  calculated  by 
o  -*  J 

N        _2 
I      (x   -  x) 

v2  =  X  "  1   \ (179) 

n  -  1 

and  for  t  distribution  with  (n  -  1)  degrees  of  freedom,  the  values  for  ±  95% 

confidence  limits,  0,  were  calculated  from 

0  =  t  v(n)1/2  .  (180) 

Values  of  t  were  taken  from  tables  in  reference  (45)  at  0.025  probability 
level  for  the  appropriate  number  of  degrees  of  freedom. 

The  statistical  curve  for  Ni  in  0.01  N  HC1-H  0  solution  obtained  in  this 
work  is  shown  in  Fig.  8  together  with  the  corresponding  curve  obtained  by 
Bockris  (43)  under  identical  conditions.   The  curves  are  in  fair  agreement 
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Table  V.   Typical  run  on  Nickel  in  0.01  N  HC1-H20 
at  25°C,  pre-elec trolysis  for  15  hours 


Current  Density 

Potential  Difference 

Overpotential 

i/A 

AV 

n 

,      2 
amps /cm 

volts 

volts 

0 

+0.141 

+0.023 

A.  80  x 

io-7 

+0.0830 

-0.035 

5.70  x 

ID"7 

+0.0770 

-0.04X 

9.44  x 

lO"7 

+0.0600 

-0.058 

1.37  x 

lO"6 

+0.0460 

-0.072 

1.94  x 

lO"6 

+0.0300 

-0.088 

2.92  x 

10"6 

+0.0140 

-0.104 

4.57  x 

lO"6 

-0.0040 

-0.122 

6.29  x 

lO"6 

-0.0260 

-0.144 

9.15  x 

lO"6 

-0.0370 

-0.155 

1.26  x 

lO"5 

-0.0520 

-0.170 

1.77  x 

lO"5 

-0.0670 

-0.185 

2.52  x 

10"5 

-0.0790 

-0.195 

3.77  x 

lO"5 

-0.0950 

-0.213 

5.14  x 

lO"5 

-0.1130 

-0.231 

6.86  x 

lO"5 

-0.1320 

-0.250 
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Table  VI,   Typical  run  on   Nickel  in  0.01  N  DC1-D  0 
at  25°C,pre-electrolysis  for  26  hours 


Current  Oensity 
i/A 

amps /cm 


Potential  Difference 

AV 


8.20  x  10 
1.00  x  10 
1.55  x  10" 
2.32  x  10" 
3.30  x  10 
5.00  x  10" 
7.80  x  10" 
1.15  x  10" 
1.52  x  10" 
2.30  x  10 
3.30  x  10 
4.40  x  10" 
6.60  x  10 
9.00  x  10 
1.15  x  10" 
1.52  x  10 
2.25  x  10" 


-6 


-6 


-5 


-5 


-o 


-D 


-4 


volts 


+0.2366 
+0.0314 
+0.0212 
-0.0011 
-0.0238 
-0.0336 
-0.0609 
-0.0769 
-0.1009 
-0.1106 
-0.1365 
-0.1455 
-0.1681 
-0.1809 
-0.1968 
-0.2138 
-0.2286 
-0.2418 


Overpotential 
n 

volts 


+0.0952 
-0.1100 
-0.1202 
-0.1425 
-0.1652 
-0.1750 
-0.2023 
-0.2183 
-0.2423 
-0.2521 
-0.2779 
-0.2869 
-0.3105 
-0.3223 
-0.3382 
-0.3552 
-0.3700 
-0.3832 
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with  each  other.   Fig.  9  shows  the  statistical  Tafel  lines  in  0.01  N  HC1-H  0 
and  0.01  N  DC1-D  0  solutions  at  25°C.   The  shaded  area  around  the  statistical 
line  is  determined  by  taking  the  standard  deviation,  a,  given  by 

n        -  2 
£   (x    x) 

o  =  ±[n=  X     i ]1/2  (181) 

n 

where  x.'s  are  the  overpotential  values  for  a  fixed  current  density  for  all 

the  Tafel  lines  under  consideration  and  x  is  the  mean  of  these  values. 

Results  obtained  on  stainless  steel  were  treated  statistically  in  the  same 

way. 


71 


Hi 

id   • 

IO 

oo 

ro 

i 

Zo 

m 

II.01 

o 

1- 

z 

O 

O^ 

i 

t 

<   d 

NO 

E 

Q. 

K) 

o 

* 

1 

oo 

ST 

o 

E 

O  z 

o 

X 

M 

r-  _ 

<o 

a. 

o  Z: 

£ 

o 

q 

a 

a: 

U) 

o 

09 

1 

— 

U-2. 

o 

< 

o» 

o 

w 

""^ 

20 

-1  X 

ro 

c 

tf> 

_l  — 

i 

UJ  w 

u. 

X 


< 

o 

r- 

co 

I 

(A 


0> 

ill 


(SHOA)^ 


72 


4.2  Measurements  on  Stainless  Steel  Electrodes 
4.2.1   Tafel  Slopes  and  Stoichiometric  Number,  u 

Typical  measurements  of  overpotential  at  various  cathodic  current  densi- 
ties made  on  stainless  steel  electrode  in  HC1  and  DC1  solutions  are  presented 
in  Tables  VE  and  VIH.   Figs.  10  and  12  show  plots  of  typical  runs  in  HC1  and 
DC1  solutions  on  stainless  steel  with  and  without  pre-electrolysis .   There 
is  a  good  amount  of  departure  from  linearity  of  Tafel  curves  for  hydrogen  and 
deuterium  runs  in  case  of  insufficient  solution  purity.   This  showed  clearly 
the  remarkable  effect  of  pre-electrolysis  on  the  data  obtained.   The  Tafel 
curves  in  such  impure  solutions  showed  two  slopes,  one  at  lower  current 
densities  and  the  other  at  higher  current  densities.   The  same  phenomena  has 
been  observed  by  Bockris  (7)  and  others  in  presence  of  impurities.   Even  with 
pre-electrolysis,  a  small  time  variation  in  the  electrode  potential  was  observed 
and  the  most  reproducible  results  were  obtained  when  the  Tafel  line  was  estab- 
lished by  the  rapid  technique  (7),  i.e.,  by  taking  the  measurements  with 
minimum  possible  delay  as  soon  as  steady  electrode  potential  was  indicated  by 
the  null  galvanometer. 

In  all  the  Tafel  lines  two  non-linear  regions,  one  at  low  current  densi- 
ties and  another  at  high  current  densities,  existed.   Fig.  13  shows  the 
statistical  Tafel  lines  for  stainless  steel  electrodes  in  HC1-H„0  solutions  of 
different  pH  values  between  0.2  and  3.8.   Table  IX  gives  the  corresponding 

values  of  b  and  log  i  .   The  ±   95%  confidence  limits  on  these  values  were 

o 

calculated  by  the  method  mentioned  above.   As  seen  from  Table  IX  no  dependence 
of  the  Tafel  slopes  on  pH  was  observed  indicating  that  the  reaction  mechanism 
did  not  change  with  the  change  in  pH.   The  variation  of  log  i  and  n  with  pH 
is  shown  in  Figs.  14  and  15. 
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Table  VII.   Typical  run  on  stainless  steel-304  in  0.01  N  HC1-H20 
solution  at  25°C,  pre-electrolysis  for  15  hours 


Current  density 
i/A 

,   2 
amps /cm 


Potential  Difference 

AV 

volts 


Overpotential 

n 

volts 


0 


8.30  x  10 


-7 


1.00  x  10 


1.65  x  10 


-6 


2.40  x  10 


-6 


3.40  x  10 


-6 


5.20  x  10 


8.05  x  10 


-6 


1.15  x  10 


-5 


1.65  x  10 


-5 


2.30  x  10 


-5 


3.50  x  10 


-5 


4.90  x  10 


-5 


7.85  x  10 


-5 


1.10  x  10 


-4 


1.75  x  10 


-4 


2.30  x  10 


3.60  x  10 


4.90  x  10 


-4 


+0.1974 
+0.0035 
-0.0020 
-0.0260 
-0.0685 
-0.0680 
-0.1020 
-0.1176 
-0.1370 
-0.1592 
-0.1780 
-0.2023 
-0.2208 
-0.2482 
-0.2670 
-0.2940 
-0.3106 
-0.3470 
-0.3540 


+0.0584 
-0.1145 
-0.1200 
-0.1440 
-0.1665 
-0.1860 
-0.2200 
-0.2356 
-0.2550 
-0.2772 
-0.2960 
-0.3203 
-0.3388 
-0.3662 
-0.3850 
-0.4120 
-0.4286 
-0.4650 
-0.4720 
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Table  VIII.   Typical  run  on  stainless  steel-304  in  0.01  N  DC1-D  0 
at  25°C,  pre-electrolysis  for  26  hours 


Potential  Difference  Overpotential 

av  n 

volts volts 

+0.2085  +0.067 

-0.0565  -0.198 

-0.0585  -0.200 

-0.0765  -0.218 

-0.0975  -0.239 

-0.1185  -0.260 

-0.1455  -0.287 

-0.1715  -0.313 

-0.1915  -0.333 

-0.2185  -0.360 

-0.2385  -0.380 

-0.2655  -0.407 

-0.2865  -0.428 

-0.3085  -0.450 

-0.3385  -0.480 

-0.3615  -0.503 


Current  Density 

i/A 

2 

amps/cm 

0 

8.10  x 

10"7 

1.00  x 

io-" 

1.65  x 

ID"" 

2.40  x 

IO"" 

3.40  x 

IO"" 

5.10  x 

IO"" 

7.95  x 

IO"" 

1.10  x 

IO"5 

1.65  x 

IO"5 

2.30  x 

io"5 

3.40  x 

io"5 

4.80  x 

io-5 

7.65  x 

io"5 

1.10  x 

io'4 

1.60  x 

IO"4 
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Fig.   14         VARIATION      OF     log|0i0      WITH      pH 
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Fig.    15         VARIATION      OF      OVERPOTENTIAL  (at   log,0lc  * -4.8) 
WITH      pH. 
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Table  X  gives  the  values  of  b  and  i   for  stainless  steel  in  H„0  and  D„0 

o  2       2 

solutions.   It  is  obvious  that  b  is  almost  the  same  in  the  two  cases,  indica- 
ting that  there  is  no  change  in  the  reaction  mechanism  due  to  deuterium 
substitution.   It  is  also  observed  that  these  values  of  b  and  those  recorded 
in  Table  IX  were  always  above  0.051.   Comparison  of  experimental  values  of  b 
with  the  corresponding  theoretically  calculated  values  given  in  Table  XI  (44) 
shows  that  the  rate  determining  step  in  the  reaction  mechanism  on  the  stainless 
steel  surface  is  either  simple  discharge  or  electrochemical  desorption.   It  has 
been  established  that  the  rate  determining  mechanism  on  Fe  is  simple  discharge 
(44)  and  that  on  Ni  is  electrochemical  desorption  (43)  both  of  which  have 
theoretical  b  values  of  0.051  (Table  XI).   Assuming  that  each  of  the  three 
elements,  Fe,  Ni,  and  Cr  composing  the  stainless  steel  alloy  will  keep  its 
individuality  as  indicated  in  Section  2.5,  it  can  be  concluded  that  the  rate 
determining  step  in  hydrogen/deuterium  evolution  on  Cr  is  either  simple  dis- 
charge or  electrochemical  desorption.   To  the  best  of  author's  knowledge  no 
direct  study  of  hydrogen  evolution  on  Cr  electrode  surfaces  has  been  made 
because  of  its  high  reactivity. 

Table  XII  includes  values  of  the  stoichiometric  number,  u,  defined  by 
Eq.  (107) 

»=4rWdic>n+o  <107) 

(dn/d  i  )      ,  the  slope  of  the  non-linear  part  of  n  -  i  curve,  is  rather 

c  n  -*■  o  c 

difficult  to  measure,  since  it  gradually  varies  and  could  give  rise  to 
indefinite  values  of  u'.   Generally  speaking  it  is  very  difficult  to  measure 
the  value  of  u  in  acid  solutions  due  to  slight  dissolution  of  the  electrode 
which  occurs  at  low  current  densities.   As  seen  in  Table  XII,  values  of  u 
vary  between  0.84  and  2.12  with  an  average  value  of  u  =  1.51.   This  is  also 


Table  X.   Statistical  values  of  i  and  b  in  HC1-H  0 
and  DC1-D  0  solutions  at  25°C 
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Solution    Normality    Temperature 


HC1-H  0 


DC1-D20 


0.01 
0.01 


-7 


25°C 

1.1019  x  10 

0.06103 

25°C 

4.3484  x  10~8 

0.06377 

2.54 


Table  XI.   Theoretical  values  of  u  and  b  for  various 
mechanisms  of  cathodic  hydrogen  evolution  in  acid  solutions 


Mechanism 


RT 


Numerical 
Value  of  b 


Slow  Discharge 

Fast  Atomic 
Desorption  (A) 


2RT/F 


0.0512 


Slow  Atomic 
Desorption 

Fast  Discharge  (B) 


RT/2F 


0.0128 


Slow  Discharge 

Fast  Electrochemical  (C) 


2RT/F 


0.0512 


Slow  Electrochemical 
Fast  Discharge  (D) 


2RT/3F 
2RT/F 


0.0171 
0.0512 
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Table  XII.   Typical  values  of  the  stoichiometric  number,  m , 
in  HC1-H  0  solutions  of  various  pH  at  25°C 


2 
PH  i   amps/cm  non-linear 

o  . 

slope  at  n^o 


0.2 

4.286 

X 

io-6 

0.2 

5.985 

X 

io-6 

0.2 

5.664 

X 

IO"6 

1.0 

1.190 

X 

IO"6 

1.0 

0.967 

X 

io"6 

1.0 

0.905 

X 

io"6 

1.5 

0.312 

X 

io"6 

1.5 

0.554 

X 

io"6 

1.5 

0.384 

X 

io"6 

2.0 

0.759 

X 

10"' 

2.0 

1.507 

X 

io"' 

2.0 

0.954 

X 

10"' 

2.5 

0.536 

X 

io"' 

2.5 

0.657 

X 

io"' 

2.5 

0.383 

X 

io"' 

5.06  x  103  1.69 

3.24  x  103  1.51 

4.81  x  IO3  2.12 


1.57  x  IO4  1.45 

2.12  x  IO4  1.60 

2.59  x  10  1.82 


5.90  x  10  1.43 

3.51  x  IO4  1.51 

5.70  x  IO4  1.70 


2.23  x  IO5  1.32 

1.39  x  IO5  1.63 

1.13  x  IO5  0.84 


2.25  x  IO5  0.94 

2.89  x  IO5  1.48 

5.34  x  IO5  1.59 
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the  most  probable  value.   Comparison  of  these  values  with  the  theoretically 
calculated  values  of  u  shown  in  Table  XII  does  not  lead  to  a  definite  con- 
clusion regarding  the  rate  determining  step  in  the  reaction  mechanism.   It 
does,  however,  favor  the  above  conclusions  arrived  at  by  comparing  the  ex- 
perimental and  theoretical  values  of  the  Tafel  slopes. 

4.2.2   Isotopic  Ratio  and  Surface  Coverage 

Fig.  16  shows  the  statistical  Tafel  lines  of  stainless  steel  in  0.01  N 
HC1-H  0  and  0.01  N  DC1-D-0  solutions.   Values  of  the  exchange  current  densi- 
ties and  the  isotopic  ratio,  R,  defined  as 

R  =  (i  )    /(i  )  (182) 

°  H20   °  D20 

are  given  in  Table  X.   R  has  been  found  to  be  characteristic  of  the  rate 
determining  step  in  the  mechanism  of  the  hydrogen/deuterium  evolution  reaction 
(44). 

Surface  coverage,  defined  as  the  fraction  of  the  surface  covered  with 
hydrogen/deuterium  atoms  may  affect  the  isotopic  ratio,  R.   Surface  coverage 
effects  can  arise  from:   (a)  differences  in  the  steady-state  surface  concen- 
trations of  adsorbed  H  or  D  in  the  rate  determining  step,  and  (b)  the 
dependence  of  heats  of  adsorption  of  H  and  D  on  coverage.   It  has  been  shown 
by  Conway  (18)  that  effects  of  type  (a)  are  small,  particularly  for  highly 
covered  surfaces.   Experimentally  (46)  at  smooth  platinum,  surface  coverage 
by  adsorbed  D  is  about  15  to  20%  greater,  at  a  given  potential,  than  that  by 
H.   In  this  case  the  rate  of  the  electrochemical  and  atom  desorption  reactions 
would  tend  to  be  slightly  greater  for  D  entities  than  for  H  entities  were  it 
not  for  the  large  opposite  effect  of  zero  point  energy  differences.   Effects 
of  type  (b)  can  be  much  larger  than  those  of  type  (a) .   Conway  (44)  found 
that  if  simple  discharge  were  rate  determining  the  value  of  R  would  be 
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independent  of  surface  coverage  while  in  case  of  electrochemical  desorption 
surface  coverage  will  have  a  direct  effect  on  R.   On  the  basis  of  the  values 
of  R  recorded  in  Table  XIII  and  on  the  assumption  that  the  overall  value  of 
R  in  case  of  stainless  steel  is  the  sum  of  the  partial  contributions  of  Fe, 
Cr  and  Ni  and  that  alloying  does  not  change  the  rate  determining  step  of 
individual  elements,  the  following  analysis  can  be  presented. 

Considering  the  rate  determining  steps  on  Fe  and  Ni  as  simple  discharge 
(47),  and  electtrochemical  desorption  (44),  respectively,  and  that  on  Cr  as 
either  simple  discharge  or  electrochemical  desorption,  as  indicated  above 
from  the  experimental  value  of  the  Tafel  slope,  b,  the  experimental  value 
of  R  can  be  used  to  identify  the  rate  determining  step  on  Cr  by  calculating 
the  theoretical  values  of  R  for  the  two  cases:  (i)  when  rate  determining  step 
on  Cr  is  a  simple  discharge  and  (ii)  when  the  rate  determining  step  on  Cr  is 
electrochemical  desorption  and  comparing  the  values  obtained  with  the  experi- 
mental value.   Table  XIV  shows  the  values  of  R  calculated  for  the  above  two 
cases.   Comparing  the  eyperimental  value  of  R  =  2.54  determined  in  this  work 
and  recorded  in  Table  X  with  the  values  in  Table  XIV  it  can  be  concluded  that 
the  rate  determining  step  of  the  hydrogen/deuterium  evolution  reaction  on  Cr 
is  simple  discharge. 


4.2.3  Difference  in  Heats  of  Activation  and  the  Limiting 
Value  of  the  Separation  Factor,  S 


As  indicated  in  Section  2.4.2  the  separation  factor,  S,  can  be  determined 
from 

tt  (H  reactants)  [F°/tt  (H  reactants)] 
S  =  (C  /C) 


D  H  I   tt  (D  reactants)  [F°/tt.(D  reactants)] 
a  D   r 


(135) 
exp[(AH*)D  -  AH*y/RT]  . 


Table  XIII.   Isotopic  ratios,  R,  of  exchange  currents 
for  various  mechanisms 
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Surface 

Coverage, 

Mechanism 

0 

0.5 

0.75 

1 

Simple  Discharge 

2.2 

2.2 

2.2 

2.2 

Electrochemical 

9.6 

6.8 

6.0 

5.6 

Atom  Desorption 

18.0 

6.0 

3.5 

2.0 

Table  XIV.   Theoretical  values  of  R  for  stainless  steel 
for  Cases  (i)  and  (ii) 


Contribution  to  R  from  the  component 


Case 


Surface 

Coverage 

Fe 

Ni 

Cr 

R 

1.00 

1.628 

.448 

.397 

2.473 

0.75 

1.628 

.48 

.397 

2.525 

0.50 

1.628 

.544 

.397 

2.569 

0 

1.628 

.767 

.397 

2.792 

(i) 


(ii) 


1.00 
0.75 
0.50 
0 


1.628 

.448 

1.009 

3.085 

1.628 

.48 

1.08 

3.188 

1.628 

.544 

1.222 

3.394 

1.628 

.767 

1.73 

4.125 

Calculated   from  Table  XIII. 
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If  the  rate  determining  step  in  hydrogen/deuterium  evolution  is  simple  dis- 
charge and  the  measurements  made  on  pure  HC1-H„0  and  DC1-D_0  solution  of 


same  acidic  concentration 


(Wi  - l 


it  (H  reactants) 

_§ =  ± 

tt  (D  reactants) 
a 


[F°/tt_(H  reactants)] 

— H — i = 

[F°/tt  (D  reactants)] 


and  Eq.  (135)  reduces  to 

S  =  exp[(AH*  -  AH*)/RT]  (183) 

where  AI-L.  and  AH  are  the  heats  of  activation  of  hydrogen  and  deuterium 
reactions  respectively.   S  can  be  directly  found  from  the  above  equation  if 


AAH 


=  AH  -  AH^  is  evaluated  theoretically  or  determined  experimentally. 
Values  of  AH  and  AH  were  determined  experimentally  in  this  investigation  on 

D         n 

stainless  steel  at  different  temperatures  in  0.01  N  HC1-H  0  (Fig.  17)  and  0.01 
N  DC1-D„0  solution  (Fig.  18)  by  using  the  relation 

log  i   =  log  B  -  AH  /2.303  RT  (184) 

where  B  is  known  as  the  Eyring's  entropy  function  (43).   The  statistical 
values  of  i  and  b  obtained  at  different  temperatures  are  presented  in  Table 
XV.   As  seen  in  the  table,  there  is  no  appreciable  change  in  Tafel  slope  due 
to  change  in  temperature  indicating  that  the  reaction  mechanism  is  independent 
of  temperature.   Figs.  19  and  20  show  the  values  of  log  i  versus  —  obtained 
in  0.01  N  HC1-H  0  and  0.01  N  DC1-D  0  respectively.   The  points  on  the  plot 
were  subjected  to  least  squares  analysis  using  —  as  the  independent  variable. 
Table  XVI  shows  the  slopes  obtained,  values  of  AH  computed  from  the  slopes, 
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Table  XV.   Statistical  values  of  exchange  current  densities,  i  , 

and  Tafel  slopes,  b,  for  stainless  steel  at  different 

temperatures  in  0.01  N  HC1-H  0  and  0.01  N  DC1-D  0  solutions 


Solution 


0.01  N  HC1-H  0 


0.01  N  DC1-D  0 


Temperature,  °C 


i  ,  amps /cm' 


i  ,  amps/cm' 


15 
25 
35 
45 
55 


0.8500  x  10 
-7 


-7 


1.101  x  10 


1.695  x  10 


-7 


2.399  x  10 


-7 


2.800  x  10 


-7 


0.06346 

0.06103 

0.06206 

0.0663 

0.07833 


2.940  x  10 

4.348  x  10 

5.545  x  10 

10.0  x  10 


-8 


-8 


0.05528 
0.06377 
0.06424 
0.06364 
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Table  XVI.   Standard  heats  of  activations  in  0.01  N  HCl-H-O 
and  0.01  N  DC1-D  0  solutions 


■k  * 

Least      AHH/D       AAHD  _  R 
Solution    Normality         B  Square 
Slope    Kcal/mole    Kcal/mole 

HCl-Ho0        0.01      2.73  x  10~3     -1.30     -5.95 


0.01      6.72  x  10  3     -1.55     -7.08 


1.13 


Table  XVII.   Isotopic  ratio,  R,  of  exchange  currents  for  Nickel 


2  (io)H3°+ 

Solution    Normality    Temperature    i  ,  amps/cm        b       R  =  -r\ — vn- n+ 

o   3 


HC1-H20        0.01  25°C       2.645  x  10  7    0.0439 

DC1-D  0        0.01  25°C       8.502  x  10~8    0.04895 


3.11 
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* 

and  the  difference  in  heats  of  activation,  AAH     .   Substituting  this  value 

D  —  H 

in  Eq.  (183)  gives  limiting  separation  factor,  S,  to  be  6.8.   Comparing  this 
value  of  S  with  the  theoretical  values  recorded  in  Table  IV  shows  that  the 
rate  determining  step  in  hydrogen/deuterium  evolution  reaction  is  a  simple 
discharge  mechanism.   It  should  be  noted  that  in  determination  of  AIL  or  AH 
using  Eq.  (184),  the  difference  in  zero  point  energy  of  the  hydrogen/deuterium 
activated  states  was  assumed  to  be  zero.   For  estimation  of  the  magnitude  of 
isotopic  differences  of  ::ero  point  energies  in  the  activated  states  for  proton 
or  deuteron  transfer,  we  can  regard  these  states  as  being  electrostatically 
intermediate  between  the  ions  (H„0  or  D„0  )  and  the  resulting  neutral  mole- 
cules (H„0  and  D?0) .   For  analogous  vibrational  modes,  the  frequencies  in  the 
ion  are  lower  than  those  in  the  corresponding  molecule  (48) .   It  is  assumed 
that  the  difference  of  zero  point  energy  of  the  non-dissociating  bonds,  be- 
tween the  initial  and  activated  states,  is  approximately  B  times  the  differ- 
ence between  the  zero  point  energies  of  initial  and  final  states  where  B, 
the  symmetry  factor,  is  equal  to  1/2.   Based  on  these  considerations  the 
relative  rates  of  hydrogen  and  deuterium  production  by  simple  discharge  can 
hence  be  lowered  by  a  factor  of  exp[630/591] ,  (44),  i.e.,  2.8  at  25°C, 

compared  with  the  values  previously  estimated  for  simple  discharge  reaction. 

6  8 
Hence  the  actual  value  of  S  =  -r^-r  =  2.43  which  is  in  excellent  agreement  with 

2  .  o 

the  experimental  value  of  R  =  2.54. 
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4.3   Conclusion 

Based  on  the  experimental  results  and  discussions  presented  above,  the 
following  conclusions  are  presented. 

1.  The  Tafel  parameters  and  the  isotopic  ratio,  R,  (Table  XVII) 
obtained  experimentally  for  Ni  in  0.01  N  HC1-H  0  and  0.01  N  DC1-D  0  solu- 
tions agree  within  experimental  error  with  the  corresponding  values  reported 
by  other  investigators. 

2.  The  average  value  of  the  Tafel  slope,  b,  for  stainless  steel  elec- 
trodes in  0.01  N  HC1-H20  solution  found  to  be  0.06103  indicated  that  the  rate 
determining  step  for  the  hydrogen  evolution  reaction  on  stainless  steel  is 
simple  discharge.   It  further  indicated  that  the  rate  determining  step  for 
the  hydrogen/deuterium  evolution  reaction  on  Cr  is  either  simple  discharge  or 
electrochemical  desorption. 

3.  Both  the  logarithm  of  the  exchange  current  density  and  the  over- 
potential  for  stainless  steel  in  HC1-H90  solutions  are  linearly  dependent  on 
the  pH,  the  slopes  of  these  lines  being  -0.871  and  -0.115  respectively. 

4.  Values  of  the  Tafel  slope  b  do  not  change  with  change  in  pH  indica- 
ting that  the  reaction  mechanism  is  independent  of  the  acidity  of  the  solution. 

5.  Values  of  the  stoichiometric  number,  y,  ranged  from  0.84  to  2.12. 
Comparison  between  these  values  and  the  corresponding  theoretical  values  did 
not  lead  to  any  definite  conclusions  as  to  the  rate  determining  step  in  the 
reaction  mechanism.   Other  investigators  (43,  47,  49)  have  been  faced  with  the 
same  difficulty. 

6.  Values  of  b  for  stainless  steel  in  0.01  N  HC1-H  0  and  0.01  N  DC1-D20 
were,  within  the  experimental  error,  the  same,  indicating  that  there  was  no 
change  in  the  rate  determining  step  of  the  reaction  mechanism  due  to  deuterium 
substitution. 
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7.  The  isotopic  ratio,  R,  was  determined  to  be  2.54.  This  value 
confirmed  that  the  overall  rate  determining  step  in  hydrogen/deuterium 
evolution  reaction  on  both  stainless  steel  and  Cr  is  simple  discharge. 

8.  The  idea  of  surface  coverage  has  been  found  useful  in  correlating 
the  experimental  and  theoretical  values  of  R. 

9.  Tafel  slopes  did  not  change  with  the  change  in  temperature  both  in 
HCl-H-O  and  DC1-D-0  solutions  for  stainless  steel  indicating  that  the  rate 
determining  step  is  independent  of  temperature. 

10.   AAH      determined  from  overpotential  measurements  in  HC1-H00  and 
DC1-D-0  at  different  temperatures  indicated  simple  discharge  as  the  rate 
determining  step  in  the  overall  reaction  mechanism  on  stainless  steel. 
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5.0   FURTHER  INVESTIGATION 

Results  obtained  in  this  work  if  combined  with  anodic  overpotential 
studies  would  lead  (as  indicated  in  Section  2.5)  to  the  determination  of 
corrosion  rates  of  stainless  steel,  thus  avoiding  the  tedious  and  time 
consuming  conventional  way  of  determination  of  rates  of  corrosion. 

Comparison  of  the  results  obtained  in  this  study  with  those  of  identical 
studies  on  irradiated  electrodes  would  help  reveal  the  effects  of  radiation 
on  corrosion  mechanisms. 

Similar  studies  on  other  alloys  and  on  their  individual  constituents 
would  help  in  clarifying  the  effect  of  alloying  on  the  electrochemical  be- 
havior of  elements. 


100 


6 . 0  ACKNOWLEDGEMENTS 

The  author  wishes  to  express  his  deep  gratitude  to  Dr,  S.  Z.  Mikhail 
under  whose  direction  this  study  has  been  accomplished.   Sincere  apprecia- 
tion is  extended  to  Dr.  W,  R.  Kimel,  Head,  Department  of  Nuclear  Engineering 
at  Kansas  State  University,  for  his  advice  and  encouragement.   Thanks  are 
also  extended  to  the  Department  of  Nuclear  Engineering  of  the  Kansas  State 
University  Engineering  Experiment  Station  and  to  the  United  States  Agency 
for  International  Development  "AID"  for  their  financial  support. 


101 


7.0   LITERATURE  CITED 


1.  U.  S.  Atomic  Energy  Commission 

Atomic  Energy  Facts 

U.  S.  Government  Printing  Office  1957 

2.  Henry,  H.  H. 

Materials  for  Nuclear  Power  Reactors 

Reinhold  Publishing  Corporation,  New  York  1955 

3.  Breden,  C.  R. 

Water  Chemistry  and  Corrosion,  Vol.  II 
Argonne  National  Laboratory  (Feb.  1963) 

4.  Uhlig,  H.  H. 

Corrosion  and  Corrosion  Control 

John  Wiley  &  Sons,  Inc.,  New  York  (March  1964) 

5.  Stern,  M. 

The  Electrochemical  Behavior,  Including  Hydrogen  Overvoltage,  of 

Iron  in  Acid  Environments 

J.  Elec.  Chem.  Soc,  Nov.  1955,  102;  609-616 

6.  Stern,  M. 

The  Relation  Between  Pitting  Corrosion  and  the  Ferrous-Ferric 
Oxidation  Reduction  Kinetics  on  Passive  Surfaces 
J.  Elec.  Chem.  Soc,  Oct.  1957,  104;  600-606 

7.  Bockris,  J.,  Pentland,  N.,  Sheldon,  E. 

Hydrogen  Evolution  Reaction  on  Copper,  Gold,  Molybdenum,  Palladium, 

Rhodium,  and  Iron 

J.  Elec.  Chem.  Soc,  March  1957,  104;  182-194 

8.  James,  F.  B. 

Anodic  Passivation  of  Steel  in  100  Percent  Sulfuric  Acid 

Corrosion  -  National  Association  of  Corrosion  Engineers,  July  1963, 

19;  238t-241t 

9.  Edeleanu,  C. 

A  Potentiostat  Technique  for  Studying  the  Acid  Resistance  of  Alloy 

Steels 

J.  Iron  Steel  Inst.,  Feb.  1958;  122-132 

10.  Stern,  M. 

The  Mechanism  of  Passivating-Type  Inhibitors 
J.  Elec.  Chem.  Soc,  Nov.  1958,  105;  638-647 

11.  Bockris,  J.,  Potter,  E.  C. 

The  Mechanism  of  the  Cathodic  Hydrogen  Evolution  Reaction 
J.  Elec.  Chem.  Soc,  April  1952,  99;  169-185 


102 


12.  Eyring,  H.,  Glasstone,  S.,  Laidler,  K.  J. 

Application  of  the  Theory  of  Absolute  Reaction  Rates  to  Overvoltage 
J.  Chem.  Phys.,  Nov.  1939,  7;  1053-1065 

13.  Glasstone,  S. 

Introduction  to  Electrochemistry 

D.  Van  Nostrand  Co,  Inc.,  New  York,  July  1947 

14.  Glasstone,  S. 

The  Elements  of  Physical  Chemistry 

D.  Van  Nostrand  Co,,  Inc.,  New  York,  Feb.  1961 

15.  Moelwyn-Hughes,  E.  A. 

Physical  Chemistry 

Pergamon  Press,  New  York,  Jan.  1957 

16.  Glasstone,  S. 

Theoretical  Chemistry 

D.  Van  Nostrand  Co.,  Inc.,  New  York,  Feb.  1961 

17.  Sherman,  J. 

Chem.  Rev.,  1932,  2;  93 

18.  Conway,  B.  E. 

The  Electrolytic  Hydrogen-Deuterium  Separation  Factor  and  Reaction 

Mechanism 

Proc.  Roy.  Soc.  (London),  1958,  A247 

19.  Dennison,  D.  M. 

Rev.  Mod.  Phys.,  1940,  12;  175 

20.  Born,  Z. 

Physik,  1920,  1;  45 

21.  Lange,  E. 

A.  Elektrochem,  1938,  44;  43 

22.  Owen,  B.  B.,  Sweeton,  F.  H. 

J.  Am.  Chem.  Soc,  1933,  55;  3504 

23.  Lewis,  G.  N.,  Doody,  F.  G. 

J.  Amer.  Chem.  Soc,  1933,  55;  3504 

24.  Harned,  H.  S.,  Dreley,  E.  C. 

J.  Am.  Chem.  Soc,  1939,  61;  3113 

25.  Krishenbaum,  J. 

Physical  Properties  and  Analysis  of  Heavy  Water 
McCraw-Hill  Book  Co.,  New  York,  1951 

26.  Gaydon,  A.  G. 

Dissociation  Energies  and  Spectra  of  Diatomic  Molecules 
London:   Chapman  and  Hall,  1953 


103 


27.  Reitz,  0.,  Forster,  E. 

Z.  Elektrochem.,  1938,  44;  45 

28.  Parsons,  R.,  Bockris,  J.  0. 

Trans.  Faraday  Soc,  1951,  47;  914 

29.  Conway,  B.  E.,  Bockris,  J.  0. 

Modern  Aspects  of  Electrochemistry,  Chap.  II 
Academic  Press,  New  York 

30.  Laidler,  K.  J. 

J.  Chem.  Phys.,  1954,  22;  1740 

31.  Papee,  H.,  Canady,  W.,  Laidler,  K.  J. 

Can.  J.  Chem.,  1956,  34;  1677 

32.  Conway,  B.  E.,  Bockris,  J.  0. 

J.  Chem.  Phys.,  1957,  26;  532 

33.  Herzberg,  G. 

Molecular  Spectra  and  Molecular  Structure 
Van  Nostrand,  New  York,  Aug.  1959,  1;  501 

34.  Bigeleisen,  J. 

J.  Phys.  Chem.,  1952,  56;  823 

35.  Bockris,  J. 

Modern  Aspects  of  Electrochemistry 

Butterworths  Scientific  Publications  (London) ,  April  1954 

36.  Topley,  B.,  Eyring,  H. 

J.  Chem.  Phys.,  1934,  2;  219 

37.  Eucken,  A.,  Bratzler,  K. 

Z.  Phys.  Chem.,  1935,  174A;  273 

38.  Okamoto,  G.,  Horiuti,  J.,  Hirota,  K. 

Sci.  Pap.  Inst.  Phys.  Chem.  Res.  Tokyo,  1936,  29;  223 

39.  Stern,  M. ,  Geary,  A.  L. 

Electrochemical  Polarization 

J.  Electrochem.  Soc,  1957,  104;  56 

40.  Mueller,  W.  A. 

Derivation  of  Anodic  Dissolution  Curve  of  Alloys  from  Those  of  Metallic 

Components 

Corrosion  Feb.  1962,  Vol.  18,  73t-79t 

41.  Uhlig,  H.  H. 

Z.  Electrochem.,  1958,  62;  700 

42.  Evans,  U.  R. 

Z.  Electrochem.,  1952,  62;  619 


104 


43.  Bockris,  J.,  Potter,  E.  C. 

The  Mechanism  of  Hydrogen  Evolution  at  Nickel  Cathodes  in  Aqueous 

Solutions 

J.  Chem.  Phys.,  April  1952,  20;  614-628 

44.  Conway,  B.  E. 

Kinetics  of  Electrolytic  H/D  Evolution 
Proc.  Roy.  Soc,  1960,  A256;  128-144 

45.  Freund,  E. 

Mathematical  Statistics 
Prentice-Hall,  Inc.  1962 

46.  Rozenthal,  K.  and  others 

Zh.  Fiz.  Khim.,  USSR,  1945,  19;  601 

47.  Bockris,  J.,  Koch,  D.  F.  A. 

Comparative  Rates  of  the  Electrolytic  Evolution  of  H„  and  D„  on  Fe,  W,  Pt 
J.  Chem.  Phys.,  Nov.  1961,  65;  1941-1948 

48.  Falk,  M.,  Giguere,  P.  A. 

Can.  J.  Chem.,  1957,  35;  1195 

49.  Azzam,  A.  M.,  Bockris,  J.  M.,  Conway,  B.  E.,  Rosenberg,  H. 

Some  Aspects  of  the  Measurement  of  Hydrogen  Overpotential 
Trans.  Faraday  Soc,  1950,  46;  918-927 


105 


APPENDIX  A 


Least  Squares  Fit  for  the  Tafel  Region  and  Calculation  of  the 
Stoichiometric  Number  from  the  Non-linear 
Region  of  the  Overpotential  Curves 


The  Tafel  equation  is 

n  =  a  +  b  In  i  .  (1) 

c 

Substituting        y  =  ri  (2) 


x  =  In  i  (3) 

c 

y  =  a  +  bx  .  (4) 

The  deviation  is   E.  =  y.  -(a  +  bx)  .  (5) 

Therefore,      (E±)2  =  [y.  -(a  +bx±)]2  (6) 

and  E(E.)2  =  E  [y.  -(a  +  bxj]2  .  (7) 

i    i=l   i         i 

By  minimizing  with  respect  to  a 

6E(E.)2 

6a    =  E[-2(y.  -  (a  +  bx±) ]  =  0 

or  =  -2E[y   -  (a  +  bx±) ]  =  0  (8) 

and  with  respect  to  b 

6E(E.)2 

*   =  E[2(y,  -  (a  +  bx.)(-x.))]  =  0 
ob  l  11 

=  -2Z[(y.  -  (a  +  bXi))x.]  =  0  (9) 

Therefore,  rewriting  equations  (8)  and  (9)  and  equating  to  zero, 

-Ey .  +na  +  bEx.  =0 
}  x  1 

or  na  +  (Ex.)b  =  Ey  (10) 
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2 
Therefore,     -Ey.x.  +  aEx.  +  bEx.   =  0 
'  i  l      l      l 

or  (Zx.)a  +  (Zx.2)b  =  Ey.x..  (11) 

i        l        i  1 

Solving  equations  (10)  and  (11)  by  Cramer's  rule 

n£y.x.  -  Ex.Ey. 

b  -  — V l-ji  (12) 

nZx.   -  (£x.) 
l       l 

2 
Ex. Ey.x.  -  Ey^x. 

a  -  ^ 2~  ^   ^ 

nix.   -  (Ex.) 


IBM-1410  Program  for  Least  Squares  Fit  for  the  Tafel 
Region  and  Calculation  of  the  Stoichiometric  Number 
from  the  Non-linear  Region  of  the  Overpotential 

Curves 

M0N$$     JOB  TAFEL  LEAST  SQUARE  SLOPE 
M0N$$     COMT  5,3  PAGES,  SHIV  N  KAUL, MASTERS  THESIS 
M0N$$     ASGN  MJB.12 
M0N$$     ASGN  MG0,16 
M0N$$     MODE  GO, TEST 
M0N$$     EXEQ  FORTRAN,,,, ,,, TAFEL 
DIMENSION  AMP(40),ETA(40) ,EAMP(A0) 

1  FORMAT (F6. 2, 313) 

2  FORMAT(E8.2,F12.6) 

3  FORMAT (5H  AMP=,E14 .4 , 5H  ETA=,F12 .6, 18H  NON-LINEAR  SL0PE=, E16 .8) 

4  FORMAT(4HKPH=,F6.2,4H  A=,E14.8,18H   ERROR  PER  POINT=,E14 .8) 

5  FORMAT (1HL,20X,10HDATA  SET   ,13) 

6  F0RMAT(5H  AMP=, E14 .4 , 5H  ETA=, F12 .6, 9H  DELTA  E=,F12.8) 

7  FORMAT (11HKINTERCEPT=,E14. 8, 28H  LINEAR  LEAST  SQUARE  SLOPE=, E18 .8) 
100  READ(1,1)  PH,NP1,NP,NSET 

IF(NP.EQ.O)  CALL  EXIT 
C      NP1  IS  NUMBER  OF  POINTS  IN  NON-LINEAR  REGION 
C      NP  IS  TOTAL  NUMBER  OF  POINTS 

WRITE  (3,5)NSET 

READ  (1,2)  (AMP(I),ETA(I),I=1,NP) 
C     CONVERT  INPUT  ETA  TO  NEGATIVE  VALUE 

DO  70  1=1, NP 
70  ETA(I)=-ETA(I) 
C      CALCULATE  DELTA  E 
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DO  50  1=1, NP 
DELE=FTA(I)+.059*PH 

50  WRITE  (3,6)  AMP (I) ,ETA(I) ,DELE 
C     CALCULATE  SLOPE  IN  NON-LINEAR  REGION 

DO  20  1=1, NP1 

SL0PE=(ETA(I+1)-ETA(I))/(AMP(I+1)-AMP(I)) 
20  WRITE  (3,3)  AMP(I) , ETA(I) , SLOPE 
C      E  PREFIX  DENOTES  QUANTITY  IN  LOGRITHMIC  DOMAIN 

NP11=NP1+1 

DO  30  I=NP11,NP 
30  EAMP(I)=ALOG(AMP(I)) 

NPL=NP-NP1 
C     NPL  IS  NUMBER  OF  POINTS  IN  LINEAR  REGION 
C     CALL  SUBROUTINE  TO  PERFORM  LEAST  SQUARE  ANALYSIS 

CALL  LSTSQ1 (ETA, EAMP , NP11 ,NP , EAO , ESLOPE) 

SL0P=1 . /ESLOPE 

A=SLOP*EAO 

AO=EXP(EAO) 
C     CALCULATE  TOTAL  POSITIVE  ERROR  PER  POINT  IN  LINEAR  REGION 

ERR=0 . 

DO  40  I=NP11,NP 
40  ERR=ERR+ABS (ESLOPE*ETA ( I ) +EAO-EAMP ( I ) ) 

ANPL=NPL 

ERR=ERR/ANPL 

WRITE  (3,4)  PH,A,ERR 

WRITE  (3,7)AO,SLOP 

GO  TO  100 

END 
MON$$      EXEQ  FORTRAN, ,,,,, ,LSTSQ1 

SUBROUTINE  LSTSQ1(X, Y,K1,K2 ,AO, Al) 

DIMENSION  X(40),Y(40) 

C12=0. 

C1K=0. 

C22=0. 

C2K=0. 

DO  10  I=K1,K2 

C12=C12+X(I) 

C1K=C1K+Y(I) 

C22=C22+X(I)**2 
10  C2K=C2K+X(I)*Y(I) 

C21=C12 

C11=K2-K1+1 

D=C11*C22-C12*C21 

A0=(C1K*C22-C2K*C12)/D 

A1=(C2K*C11-C1K*C21)/D 

RETURN 

END 
MON$$      EXEQ  LINKLOAD 

CALL  TAFEL 
MON$$      EXEQ  TAFEL, MJB 
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APPENDIX  B 


Chemical  Analysis  of  the  Stainless  Steel  Sample  Used 
for  Preparation  of  Stainless  Steel  Electrodes 


Name  of  the  supplier:     American  Steel  and  Wire  Division  of  United  States 

Steel  Corporation 

Description  of  Material:   USS10-8  Type  304  Stainless  Steel  Wire  BRT  ANLD 

Gray  Matte  Finish 

Diameter:  1/16" 

Analysis:                Carbon  .08% 

Manganese  .94 

Phosphorus  .029 

Sulphur  .018 

Silicon  .53 

Nickel  8.76 

Chromium  18.76 
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APPENDIX  C 


Report  on  the  Spectroscopically  Pure  Nickel  Sample 
Used  for  Preparation  of  Nickel  Electrodes 


Name  of  the  supplier:      Johnson,  Matthey  &  Co.,  Limited,  London 

Description  of  Material:   Nickel  sheet,  1/16"  (The  material  has  been  prepared 

in  special  Chemicals  and  Metals  Laboratory  and  of  a 
high  degree  of  purity.) 

Spectrographic  Examination:   A  spectrographs  examination  was  made  by  means 
of  a  constant  current  D.C.  arc,  taking  5.6  amps.,  between  pure  graphite  elec- 
trodes.  Weighed  quantities  of  the  sample  were  arced  in  thin-walled  anode  cups 
against  a  horizontal  machined  electrode  as  the  cathode,  both  electrodes  being 
water-cooled. 

Spectra  were  photographed  simultaneously  on  an  Ilford  Long  Range  Spectrum 
plate  with  a  Medium  Spectrograph  and  an  Ilford  Ordinary  plate  with  a  Littrow 
Spectrograph  covering  the  wavelength  range  2200-3700A  in  two  exposures. 

Estimates  of  the  quantities  of  impurities  present  were  made  by  visual 

comparison  of  the  spectra  with  those  of  synthetic  standards,  arced  in  a  manner 

similar  to  that  used  in  the  test. 

Element  Estimate  of  Quantity  Present 

parts  per  million 

Iron  8 

Silicon  8 

Aluminum  5 

Copper  5 

Calcium  2 

Magnesium  2 

Manganese   each  element  1 

Silver     less  than 
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The  following  elements  were  specifically  sought  but  not  detected,  i.e. 
either  they  are  not  present  or  they  are  below  the  limits  of  detection  by  the 
described  examination  procedure:   As,  Au,  B,  Ba,  Bi,  Cd,  Co,  Cr,  Cs,  Ga,  Ge, 
Hf,  Hg,  In,  Ir,  K,  Li,  Mo,  Na,  Nb ,  Os,  P,  Pb,  Pd,  Pt,  Rb,  Re,  Rh,  Ru,  Sb,  Se, 
Sn,  Sr,  Ta,  Te,  Ti,  Tl,  V,  W,  Zn,  Zr. 

During  fabrication  some  surface  contamination  is  unavoidable  and  for  this 
reason  the  samples  used  for  this  spectrographs  analysis  have  been  pickled  in 
hot  dilute  nitric  acid  before  examination. 
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APPENDIX  D 


Calculation  of  the  Difference  in  the  Isotopic  Heats  of  Activation 
from  Potential  Energy  Diagrams 


The  initial  and  final  state  potential  energy  curves  for  H  and  D  are  shown 
in  figure  (21)  : 


INTERATOMIC   DISTANCE 


Fig.  (21)   Potential  Energy  Diagram 

Neglecting  the  zero  point  energy  differences  of  the  initial  state,  the 
following  geometric  relationships  can  be  obtained-  from  the  above  figure. 


AH   =  a-b-BQ. 
H         1 


AHD  =  a-b+3Qf-Qf 

=  a-b-(l-3)Qf 

where  3  is  the  symmetry  factor  and  Q  and  Q  are  the  differences,  between  the 
minima  of  the  potential  energy  curves  of  hydrogen  and  deuterium  species  for 
the  initial  and  final  states,  respectively. 
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Thus  AHD  -  At^  =  3Q.  -  (1  -  B)Qf 

and  if  6  =  1/2,  AH*  -  AH^  =  l/2(Qi  -  Qf)  =  1/2  AQ 

Now, correcting  for  the  zero  point  energies  of  the  initial  state,  the  required 
activation  energies  are 

*  * 


AH 
o 


,H  '  AHH  "  <z'P-e'>H 


AHo,D  =  %  ~  (2-P'e-)D 

or  AH*   -  AH*   =  1/2  AQ  +  (z.p.e.)H  -  (z.p.e.) 

0,L)o,H  n  JJ 

where  z.p.e.  =  zero  point  energy. 
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ABSTRACT 

The  purpose  of  this  work  was  to  study  the  electrode  kinetics  of  hydrogen 
evolution  on  stainless  steel-304,  to  determine  the  effect  of  deuterium 
substitution  on  the  reaction  mechanism  and  to  evaluate  the  hydrogen-deuterium 
separation  factor  on  stainless  steel-304,  a  factor  of  extreme  importance  in 
the  electrolytic  production  of  heavy  water.   The  reaction  mechanisms  were 
studied  by  conducting  cathodic  polarization  runs  on  stainless  steel-304 
electrodes  in  ultra-purified  systems  of  0.01  N  HC1-H  0  and  0.01  N  DC1-D20. 
Comparison  between  the  experimental  values  of  Tafel  parameters  and  their 
corresponding  theoretical  values  led  to  determination  of  the  rate  determining 
step  of  hydrogen/deuterium  evolution  reactions.   The  separation  factor,  S,  was 
evaluated  both  by  determining  the  isotopic  ratio,  R,  of  exchange  current  densi- 
ties in  light  and  heavy  water  and  by  determining  the  isotopic  difference  in 
heats  of  activation  from  overpotential  measurements  carried  out  at  different 
temperatures  in  light  and  heavy  water.   The  value  of  S  evaluated  from  the 
isotopic  ratio  was  2.5  ±0.9  and  the  value  of  S  based  on  heats  of  activation 
was  2.4  ±  an  estimated  standard  deviation  of  0.9. 

Cathodic  polarization  measurements  were  also  taken  by  using  nickel 
electrodes  in  the  same  solutions.   The  Tafel  parameters  and  isotopic  ratio,  R, 
agreed  with  those  reported  by  other  investigators. 

The  values  of  Tafel  slopes  obtained  in  the  light  and  heavy  water  systems 
were  in  agreement  within  experimental  error  and  both  values  characterized  the 
overall  rate  determining  step  in  the  hydrogen/deuterium  evolution  reaction  as 
simple  discharge.   The  stoichiometric  number,  u,  was  determined  from  the  slopes 
of  the  non-linear  region  of  the  Tafel  curve  in  the  light  water  system. 


The  isotopic  ratio,  R,  was  found  to  be  equal  to  2.5  ±  0.9.   Comparison 
of  this  value  with  the  corresponding  values  obtained  from  thermodynamic 
considerations  confirmed  that  the  rate  determining  step  in  hydrogen/deuterium 
evolution  reaction  is  simple  discharge. 

Based  on  the  above  values  of  S  and  R  it  was  shown  that  the  rate  deter- 
mining step  of  the  hydrogen/deuterium  evolution  reaction  on  chromium  was  also 
that  of  simple  discharge.   This  method  made  it  possible  to  study  such  elements 
as  chromium  which  otherwise  cannot  be  studied  independently  because  of  their 
high  reactivity. 


